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Abstract 

We consider top-quarks produced at large energy in e + e~ collisions, and address the question of 
what top-mass can be measured from reconstruction. The production process is characterized by 
well separated scales: the center-of-mass energy, Q, the top mass, m, the top decay width, 1^, and 
also Aqcd; scales which can be disentangled with effective theory methods. In particular we show 
how the mass measurement depends on the way in which soft radiation is treated, and that this can 
shift the mass peak by an amount of order QAqcd/w- We sum large logs for Q 3> m 3> T t > Aqcd 
and demonstrate that the renormalization group ties together the jet and soft interactions below 
the scale m. Necessary conditions for the invariant mass spectrum to be protected from large logs 
are formulated. Results for the cross-section are presented at next-to-leading order with next-to- 
leading-log (NLL) resummation, for invariant masses in the peak region and the tail region. Using 
our results we also predict the thrust distribution for massive quark jets at NLL order for large 
thrust. We demonstrate that soft radiation can be precisely controlled using data on massless jet 
production, and that in principle, a short distance mass parameter can be measured using jets with 
precision better than Aqcd- 
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I. INTRODUCTION 



The top quark is the heaviest known fermion of the Standard Model and couples strongly 
to the Higgs sector. The most recent CDF and D0 measurements obtained a top mass, 
m t = 170.9 ±1.8 GeV [l|], with ~ 1% uncertainty. For the standard model a precise top mass 
determination is important for precision electroweak observables which test the theory at 
the quantum level, and which constrain extensions of the theory such as supersymmetry. In 
Ref. Q we derived a factorization theorem for the invariant mass distribution of high energy 
top jets for e + e~ — > ti, which allows in principle a determination of m t with uncertainty 
better than Aq C d- Such accuracy is possible because the factorization theorem separates the 
perturbative and non-perturbative contributions in terms of field theory Wilson coefficients 
and matrix elements. A virtue of our approach is that the non-perturbative matrix elements 
are universal and in some cases are straightforward to extract from other processes. In 
addition the factorization theorem provides a unique prescription for determining the Wilson 
coefficients and perturbative matrix elements at any order in the a s expansion. This level 
of control allows us to make stable predictions for the invariant mass distribution in terms 
of a short- distance top quark mass, which is not limited in precision by Aqcd- 

Determining the top mass with jet reconstruction methods in general faces issues such 
as i) defining an observable that is sensitive to the top mass, ii) soft gluon interactions and 
color reconnection, iii) uncertainties from higher order perturbative corrections, iv) the large 
top quark width ~ 1.4 GeV, and other finite lifetime effects, v) final state radiation, 
vi) initial state radiation, vii) treatment of beam remnants, viii) underlying events, and ix) 
parton distributions. In Ref. we addressed the definition of a suitable top quark mass 
m and issues i) through v) in the framework of electron-positron collisions at high energies 
Q ^> m, where Q is the center of mass energy and m is the top mass. 1 The analysis is 
suitable for a future linear collider. Issues vii) through ix) are avoided by treating the e + e~ 
initial state, but are important in a hadron collider environment like the Tevatron or LHC. 
Issue vi) is also greatly simplified in e + e~ annihilation, since the inclusion of initial state 
photon radiation mainly shifts Q and thus has very little impact on our analysis. 

Our analysis of top jets uses effective theory techniques to exploit the hierarchy of scales 
Q ^> m ^> T > Aqcd; and separate dynamical fluctuations. This hierarchy provides a 
systematic power counting in m/Q and T/m, and gives a clear interpretation to elements in 
the factorization theorem. In Ref we focused on developing the formalism and describing 
the main conceptual points in the factorization theorem for the invariant mass distribution 
in the peak region. The same formalism also yields a factorization theorem for the invariant 
mass distribution in the tail region above the peak. Here we use models for the soft function 



We will use m for the top mass when it is not necessary to specify the precise scheme which defines this 
parameter. 
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that are consistent for both the peak and tail regions, and carry out detailed calculations of 
perturbative quantities in the factorization theorem. We verify that the matching conditions 
which define the Wilson coefficients at the scales Q and m are infrared safe, compute one- 
loop perturbative corrections to the matrix elements, and carry out the next-to-leading-log 
renormalization group summation of large logs. For the peak region these are logs between 
the scales Q, m, T, and Aqcd, while away from the peak they are between Q 2 , m 2 , and the 
variables M 2 — m 2 and M 2 — mf described below. 

As an observable sensitive to the top mass, we considered in Ref. the double differential 
invariant mass distribution in the peak region around the top resonance: 



d 2 a 



dM 2 dM? ' 



M 2 t - m 2 ~ m T < m 2 , (1) 



where 



M t 2 =($>f), (2) 

Here X t and Xi represent a prescription to associate final state hadronic four momenta to 
top and antitop invariant masses respectively. For simplicity we call X t ^ the top and antitop 
jets, and M t j the invariant mass of the top and antitop jets respectively. The distribution 
in Eq. (|l|) has a width r ~ T t + QAqqb / Tn which can be larger than the top quark width 
Ft- The restriction M 2 ^ — m 2 ~mT<C m 2 defines the peak region, which is the region most 
sensitive to the top quark mass m. Here the dynamics is characterized by energy deposits 
contained predominantly in two back-to-back regions of the detector with opening angles of 
order m/Q associated with the energetic jets or leptons coming from the top and antitop 
decays, plus collinear radiation. The region between the top decay jets is populated by soft 
particles, whose momentum is assigned to one of M 2 or M 2 . The tail region is defined by 
invariant masses starting just past the peak where the cross-section begins to fall off rapidly, 
namely where m 2 ^> M% — m 2 and either M 2 ^ — m 2 > mT or M 2 ^ — m 2 > mT. Farther 
out, when M 2 ^ — m 2 ~ m 2 , we have an ultra-tail region where the cross-section is very small. 
We do not consider the region where M 2 ^ ~ Qm. The observable in Eq. ([I]) in the peak and 
tail regions is the main focus of our analysis. We also briefly consider the cross-section in 
the ultra-tail region. 

The result for the double differential cross-section in the peak region to all orders in a s 
is given by @ 

da ( Q 



dM 2 dM? 



°~o Hq(Q, fi m )H m (mj, , fi m , fi 

x / di + d£~B + (it - — , r t , fi) £L (s- t - 9L r t , fj) s(£ + , e~,n) 

+ 0{ ^ ) + ^ ) + Q + 0{ ^y (3) 
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FIG. 1: Sequence of effective field theories used to compute the invariant mass distribution. 



where, as indicated, power corrections are suppressed by a s m/Q, m 2 /Q 2 , T t /m, or s t f/m 2 . 
Here mj is the short-distance top quark mass we wish to measure, and for convenience we 
have defined 

s t M 2 - m 2 _ s- t M 2 t -m 2 

s t — — , S t — — , (4 J 

vtij mj mj mj 

where s t) f ~ T are of natural size in the peak region. In Eq. (|3]) the normalization factor <jq 
is the total Born-level cross-section, the Hq and H m are perturbative coefficients describing 
hard effects at the scales Q and mj, B± are perturbative jet functions that describe the 
evolution and decay of the the top and antitop close to the mass shell, and S is a nonpertur- 
bative soft function describing the soft radiation between the jets. To sum large logs B± and 
S will be evolved to distinct renormalization scales /x, as we discuss in section [11 C] below. 
For the tail region Eq. (^j becomes 



so the only changes are that the soft-function S = S paTt (£ + , £~ , fj) becomes calculable, and 
we have an additional O(Aqcd0 I ' s t,l) nonperturbative correction from the power expansion 
of the soft-function which we will include in our analysis. The result in Eq. was derived 
by matching QCD onto the Soft Collinear Effective Theory(SCET) @, 1, | |, 0] which 



in turn was matched onto Heavy Quark Effective Theory(HQET) || || |10|, |TT|, [12], [13 
generalized for unstable particles 0, [15], [16], [T^| as illustrated in Fig. |I]. The decoupling of 
perturbative and nonperturbative effects into the B± jet functions and the S soft function 
was achieved through a factorization theorem in SCET and HQET, aspects of which are 
similar to factorization for massless event shapes fl8 , 19 , [20| , |2l| . The result in Eq. (§) is an 



event shape distribution for massive particles, and can be used to determine common event 
shapes such as thrust or jet-mass distributions. Note that a subset of our results can also 
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be used to match results with the event shape cross-sections for massless jets, namely by 
using our SCET ultratail cross-section and taking the limit m — > 0. 

In general the functions B± and S depend on exactly how M t and Mj, or equivalently 
X t and Xii are defined. The factorization theorem in Eq. (|3|) holds in the form shown when 
all the soft radiation is assigned to either X t or Xj, and the probability of radiation being 
assigned to X t or Xj increases to unity when we approach the top or antitop direction [0 . 
Finally, the definition should be inclusive in the hard jets and leptons from the top decay. 
One possibility for defining Mf^ in Eq. is a hemisphere mass definition, where X tl Xi 
contain everything to the left or right of the plane perpendicular to the thrust axis. In 



this case our S is identical to the soft function of Refs. |L8], [19], £2| that appears in the 
factorization theorem for massless event shapes in the dijet region. For studies of soft- 
functions in massless event shapes see Refs. jrj| [19], ^(J, [21], |23|, |4], |25], |27], |29" 



The B± are inclusive in the jets from the top decay and collinear radiation and can be 
defined by forward matrix elements @. Other definitions to associate all soft radiation to 
the top and antitop jets can be used which modifies the required S function, but for the 
class of masses defined above leaves B± unchanged. 

The use of a short- distance mass definition in the B± jet function and a short-distance gap 



parameter in the soft function S [|30| are crucial for obtaining predictions that remain stable 
when higher order perturbative corrections are included. In Ref. [[J we showed that suitable 
mass schemes for reconstruction measurements can only differ from the pole mass by an 
amount 5m ~ T t a s , and we proposed a jet-mass scheme which satisfies this criteria. We will 
refine the criteria for this jet-mass scheme here. In Eq. the jet-mass, mj, only appears 
in the calculable Wilson coefficients and jet functions B±. The greatest sensitivity to mj 
is in B±. Through these jet functions, mj influences the spectrum of the mass distribution 
in the peak region. The spectral distribution and location of the peak are also affected by 



nonperturbative effects in the soft function S. In Ref. |30[ a gap parameter scheme based on 
moments of the partonic soft function was devised to avoid perturbative ambiguities in the 
definition of the partonic endpoint where the variables in Eq. (JJ) approach zero. Methods 
for using Eq. (H) to extract mj are discussed in detail in Ref. @. 

In this paper we determine the functions Hq, H m , B± at one- loop order in a s , and carry 
out the summation of large logs between the scales Q > m > T in Eq. (|3|). The derivation 
of results for the top jet-mass scheme are discussed in detail. We also show that there are 
constraints on the allowed soft functions, and implement a consistent method to include 
perturbative corrections in S. In our numerical analysis we extend the work in Ref. 
to one-loop order, including the summation of the next-to-leading order logarithms using 
renormalization group (RG) evolution in effective field theories. Our analysis of the tt jet 
cross section at this order includes both invariant masses in the peak region and the tail 
region above the peak, and the final results are analytic up to integration over the soft- 
function model. 



7 



For massless jets there has been a lot of work done on the program of resumming logs in 
event shape variables |3IJ gl] JJ^ |§ g§ |§ [§, g7|, |3|, 0, 0, ^ ||, . In this paper we 



do not use the traditional approach to resummation, but rather an approach that sums the 
same large logs based on the renormalization of operators in effective field theories, including 
HQET and SCET [7|. The effective theory resummation technique has the advantage of 
being free of Landau-pole singularities P5| , fi6| , since it only depends on the evaluation of 
anomalous dimensions at perturbative scales. This technique can also be extended in a 
straightforward manner to arbitrary orders, N fc LL in the resummation |57], [53]. A recent 
application of the SCET technique is the resummation for thrust in e + e~ to massless jets 
at NLL order 



In our log-summation there is an important distinction between large logs which affect the 
overall cross section normalization, and large logs that change the shape of the distribution 
in M t 2 j. In predicting the normalization in the dijet region we must sum up a series of 
double Sudakov logarithms that occur for Q ^> in and for m ^> T. However, it turns out 
that the same is not true for logs affecting the shape of the invariant mass spectrum. As we 
discuss in detail, the form of the spectrum is protected from large logs below the scale Q 
until we reach the fundamental low energy scale governing the dynamics of either the soft or 
jet functions. This conclusion is not affected by the mass threshold at m, and is valid to all 
orders in perturbation theory (ie. for both leading and subleading series of logarithms). In 
order for this cancellation to occur it is important that the invariant mass definition includes 
soft radiation at wide angles. The hemisphere mass definition of M t and Mj, as well as other 
definitions which associate wide angle soft radiation to both X t and Xj, are in this category. 
In the effective field theory this protection against the appearance of shape changing large 
logs is described by a set of "consistency conditions". From our analysis we find that the 
only shape changing large logs occur between the low energy scale fi ~ QA/m + T t where 
logs in the jet functions are minimized, and a perturbative low energy scale /i <^ A + mT t /Q 
where logs in the soft function are minimized. Here A ~ 0.5 GeV is the hadronic scale 
where the interactions are non-perturbative. As indicated there are two scales appearing in 
each of these functions, and the question of which dominates depends on the size of these 
parameters. 

The program of this paper is as follows. In Sec. [HA] we review the formulation of the 
factorization theorem for the invariant mass cross-section from Ref. 0. In Sec. [11 B| we 
show that the finite lifetime effects can be treated as a convolution of B± jet functions for 
stable top quarks with a Breit-Wigner, and we describe models for 5* that are consistent 
in the presence of perturbative corrections. In Sec [11 C| we discuss the structure of large 
logarithms and present the factorization with log resummation. In section [ITT] we discuss the 
connection between renormalization and the resummation of large logs in SCET and HQET, 
derive the consistency conditions, and summarize results for the NLL renormalization group 
evolution. Results for the matching, running, and matrix elements in SCET including the 



soft hemisphere function are given in Sec. |V[ In Sec. |V| we give matching, running, and 



matrix element results in HQET. A short-distance jet-mass scheme is discussed in detail 
in Sec. [VT], including its relation to other schemes at one-loop. In Sec. |V11| we discuss the 



non-perturbative soft function model and the scheme for the gap parameter we use in our 
numerical analysis. An analysis of the one-loop cross-section with NLL log-summation is 
given in Sec. |V111| for both the peak and tail regions. Conclusions are given in section [IX. 
Additional computational details are given in the Appendices [A]-[F]. 

II. FORMALISM 

A. Invariant Mass Cross-Section 

In this section we review the main definitions of effective theory objects needed for our 
calculations of terms in the factorization theorem in Eq. (|3|). Further details can be found 
in Ref. ||. Starting from QCD, the two-jet cross section a(e + e~ — > 7*, Z* — > j t ji) can be 
written as 

a = f>7r) 4 5 4 ( ? -p*) £ ^ (0\j}"(0)\X)(X\jr(0M, (6) 

X i,j=a,v 

where q = p e - +p e +, and q 2 = Q 2 , and L l j v is the leptonic tensor including vector and axial 
vector contributions from photon and Z boson exchange. This result is valid to all orders in 
the QCD coupling but lowest order in the electroweak interactions. The superscript res. on 
the summation symbol denotes a restriction on the sum over final states to the kinematic 
situation given in Eq. ([I]). The QCD top quark currents are J? = ■0(x)rf^(x), where 
= 7^ and = 7^75. In Ref. || we started with Eq. (|6]) and derived the factorization 
theorem for the double differential invariant mass distribution in the peak region in Eq. ([$]). 
There the factor ctq is the tree-level Born cross-section, 



a = N, 



Ana 



2Q 2 v e v t e t , Q 4 (v 2 + a 2 e )(v 2 + a 2 ) 



+ 



(7) 



Q 2 -m| ' (Q 2 -m|) 2 

wheref/ = (T^ —2Q f sin 2 9 W ) / (2 sin 9w cos 9w) anda/ = T 3 ^/(2sin6 l vi/cos^ H /). Equation (g) 
can be easily generalized to include the angular distribution in cos(#) where 9 is the angle 
between the top jet direction and the e~ momentum: 

d 3 a ctq(^) d 2 & 



dM 2 dMl d cos(9) " ctq dM 2 dM 2 t ' 



(8) 



where 

dao irN c a 2 



dcos(9) 2Q 2 



2 2Q 2 v e v t e t Q 4 (v 2 + a 2 e )(v 2 + a 2 ) ■ 
6t Q*-m 2 z + (Q 2 ~m 2 z y <^"^'^ 



f AQ 2 e 2 t a e a t 8Q 4 a e v e a t v t } ' 
+ \Q2_ m | (Q2- m 2 z ) 2 J CO 



(9) 
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The remaining functions in Eq. ([3]) include Hq, a hard-function that encodes quark-gluon 
interactions at the production scale Q, H m which encodes perturbative effects 2 at the scale 
m, B±, the jet functions for the top jet and antitop jet respectively, and the soft function 
S which encodes non-perturbative information about soft hadrons radiated between the 
hard jets. The convolution with the soft function causes a correlation between the two-jet 
functions and affects the invariant mass spectrum. Each of the functions H, H m , B + , B-, 
and S in Eq. (|3|) can be defined as matrix elements of operators in an appropriate EFT, or 
as matching coefficients between two EFT's. At the scale Q the matching of QCD currents 
onto SCET is given by a convolution formula 

^ M (0) = JdudvC{u,u,fj,)J?{uj,u,fj,), (10) 

where C contains short-distance dynamics, while Jf describes all scales that are longer 
distance than Q. After making a field redefinition || the SCET production current at 
leading order in A is given by 

where we have collinear fields and Wilson lines defined in the jet-fields Xn,w(0) — ~ 
n-V){W^ n ){0) and Xn,d°) = ~ n ' ^)(^Cn)(0), as well as soft Y- Wilson lines and 
mass-mode S- Wilson lines to be discussed below. Here the (0) indicates that the fields are 
at coordinate x M = 0; recall that this dependence carries information about residual 
momenta at scales < QX 2 = m 2 /Q. The dependence on larger momenta is encoded in 
the labels of the collinear fields |J. For example, 5(u — n ■ P) forces the total minus-label- 
momentum of {W^n) to be uj. In terms of C defined in Eq. (|T0|) , the hard-function appearing 
in Eq. (^) is simply 

H Q (Q^) = \C(Q,-Q,fx)\\ (12) 

and after including RG-evolution we have Hq^Q^h) = Hq(Q , H-q)Uh q (Q , [J>q, A*) where Uh q 
is the evolution kernel discussed below in sections [II C| and |TTT| . 

We obtain the SCET two-jet cross section by replacing the QCD current in Eq. (|(J) with 
the SCET current. The resulting expression can be factorized as discussed in Ref. 

a H Q (Q,n)M{mj,n) (13) 



dM 2 dM 2 t 



x di + di J n {s t - Q£ + ,mj,r t ,[i)J n (si- Qt ,mj,T t ,fi)S(£ + ,£ ,n,mj) 



2 The coefficient H m is also sensitive to the ratio m/Q through its anomalous dimension. Here m/Q is the 
cusp angle by which the heavy-quarks are off the light-cone j5(J. See also section ^ below. 
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This result can be used to compute the cross-section in the ultra-tail region, where s t) i ~ m 2 - 
Due to the large suppression this region is not interesting experimentally, however we will still 
discuss formal aspects of Eq. (|13|) in detail because it is an important step towards deriving 
the peak region factorization theorem, and is also important for making the analogy with 
massless event shapes. The soft function S(£ + ,£~ , fi,mj) in Eq. (|T3"| ) is the same as the 
soft function in Eq. (|), up to perturbative effects due to top-quark vacuum polarization 
graphs denoted by the extra argument mj. It can be either derived by using eikonal Ward 
identities [ETJ or properties of the coupling of usoft gluons to collinear particles in SCET [ITS 



For the case of hemisphere invariant masses it is S(£ + ,£ ) = Sh emi (£ + ,£ ) where 
S hemi (£ + ,£~) = 7^ 5^ - - <0 1 (l^) ce (0) |^> | (T^t )^ (0) |0> 

= 7^<°l ( ^ )crf (^) ce (°)^( £+ -C^)0^( £ - J V)(^ + ) e - f (^i)^(o)|o> . (14) 

The same function Shemi appears in event shapes for massless two-jet production, and besides 
the mj and T t dependence, Eq. ([13|) is analogous to the factorization theorem for massless 
dijets ^(J In Eq. ( |i~4"l) c, d, e, / are color indices, N c = 3, and the soft Wilson lines 

are 



-x. 



Y n (x) = P exp ^ — ig J dsn-A s (ns+x)j , Y%(x) = P exp (ig J dsn- A s {ns + : 
Y n \x) = P exp (ig J dsn- A s {ns+x) s j , Y, R {x) = P exp ^ — ig J dsn- A s {ns+x)^ , (15) 

with A^ = T A A A for the antitriplet representation, where T A = —{T A ) T . In Eq. ( |14"D k a s 
is defined as the soft momentum components from the state X s that are included in the 
experimental determination of M t (and k b s for Mi). We also have operators = n- P a and 
P h ~ = n ■ Pb that project out the soft momentum components kf a and k~ b 

P a \X s ) =K\X S ), P b \X s ) = k b s \X s ). (16) 

For hemisphere masses the operator P+ is defined to project out the total plus-momentum 
of soft particles in hemisphere-a (and P b ~ the minus-momentum in hemisphere-6). In Ref. |2j 
it was shown that S(£ + ,£~) does not depend on the top quark width, and when we pass 
below the top quark mass scale is only modified by a perturbative prefactor, 

S(£ + , £~,fi, mj) = T ( mj , /i) S(£ + , £~,fj). (17) 

The matching coefficient To(mj,fi) is induced by the coupling of A^ gluons to top-vacuum 
polarization bubbles at zero- momentum. This result applies at any order in a s , but at NLL 
order T = 1. 
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In Eq. ( |T3"D the mass-mode function A4(mj,fi) contains virtual perturbative corrections 
due to gluons Af^ and quarks ip m with momenta p M ~ (mj, mj, mj), and is given by 



M{m J , f i) = ± i \{0\S?S*\0)\ 



(18) 



The definition of these mass-mode S- Wilson lines is identical to those in Eq. ([T^), except that 
they involve gluon fields which couple to massive top-quarks for any momentum, and 
which have zero-bin subtractions to avoid double counting the momentum region accounted 
for by the A^ s gluons. This implies that Ai(mj, fi) only gets contributions from graphs with 
a top-vacuum polarization bubble J52], |53], |54| coupling to the Af^ gluons. At NLL order 
the function A4(mj,fi) = 1, but is relevant at NNLL order and beyond when considering 
virtual top loops. Note that due to the invariant mass constraint s <C mQ, the ip m quarks 
never appear in the final state. This is important for the validity of Eq. ([13]). 

Matrix elements of top quark collinear fields in SCET give the jet functions J n for the 
top quark jet, and J n for the antitop jet, 



Jn(Qr^-m 2 j,mj,T t ,/j,) 



1 



4nN c Q 
1 

4nN c Q 



Im 
Im 



ijd'xe— (0|T{x„ lQ (0)te(x)}|0> 
i fd'xe— (0|T{x fli _ o (x)te(0)}|0> 



(19) 



These jet functions J n and J n depend on both the mass and width of the top quarks. The 
matrix elements of collinear fields are defined with the zero-bin subtractions [ESj, which 
avoids double counting the soft region. 

For predictions in the peak region, the J n and Jn functions should be factorized further by 
integrating out the top quark mass. This is accomplished by matching onto jet functions B± 
in HQET with boosted heavy quarks. The relevant Feynman rules are given in Appendix ||. 
The jet function matching takes the simple form || 

r 



(20) 



J n (s t} mj,T t ,fi Q ) = T + (m J} fi Q )B + (s tl T t ,^ Q ) + o(— ) + o(— 



Jn(st,mj,T t ,fi Q ) = T_(mj, fi Q )B_(sf,T t , /i Q ) 



The HQET jet functions B + and S_ also depend on the residual mass term Smj that fixes 
the mass definition in HQET. They are defined by 

B±(s,T t ,ti) = lm[B±(s,T t ,n)], (21) 

where the B± are vacuum matrix elements of T-products of HQET operators 

B + (2v + -r,T t ,fi) = Jd 4 xe^ (0\T{h v+ (0)W n (0)Wl(x)h v+ (x)}\0) , 

B_(2v_.r,r t ,/i) = 7 ^ 77r [d*xe ir - x (0\T{h v _(x)W n (x)W n '(0)h v _(0)}\0) . (22) 



4nN r m 
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Here for B + (s t ,T t , //) we have s t = 2v + ■ r, while for £>-(.§?, T t , /i) we have % = 2t>_ • r. 
The gluons in W n and W 7 ^ and HQET fields h v± are only sensitive to fluctuations below m 
and are built of gluons A± describing low energy fluctuations down to p 2 ~ T 2 in the top 
and antitop rest frames respectively. In Ref. these gluons were called ultracollinear. We 
emphasize that to make the matching consistent, the collinear gluons in Eq. ( P2"| ) have zero- 
bin subtractions for the same region as those in the SCET jet functions. These subtractions 
ensure that the B± jet-functions do not double-count the soft region encoded in S, and are 
critical for ensuring that the functions B± are IR-finite, as we discuss further in Appendix [A|. 
In Eq. (|22|) the Wilson lines are 

W^x) = P exp (ig J dsn-A + (ns+x)^j , W n (x) = P exp ^— ig J dsn-A + (ns+x)j, (23) 

with analogous formulas for W n and W\ in terms of n ■ A_. Note that if the Wilson lines 
W n and W n were absent, then B± would just define the HQET heavy quark/ ant iquark 
propagators ||. The Wilson lines, let's say for B + , encode the color dynamics of gluons that 
are soft in the top rest frame and come from the highly boosted antitop quark, and they 
render this vacuum matrix element into a gauge-invariant physical object. The analogous 
situation with top and antitop switched applies for the vacuum matrix element For the 
SCET jet functions, the Wilson lines appearing in Eqs. (|IT]) and (|T^) have the analogous 
physical interpretation where the top quark mass has not yet been integrated out. 

In the final factorization theorem in Eq. (^|) we have the B± functions, as well as the 
matching condition for the mass fluctuations, 

if m (m,/i) = T + (m, yu)T_(m, n)T (m, /i)M(m, /i) . (24) 

In bHQET all dynamic effects associated with the top-quark mass appear in H m , and 
there are no mass- mode quarks or gluons in this theory. Since To and M. (m, fi) encode 
finite matching corrections at the scale fi ~ m due to top-vacuum polarization, we have 
T (m, fi)A4(m, fi) = 1 + (D(a 2 ), and so these factors drop out from our NLL analysis. 
Therefore in later sections we simply use H m = T + T_. This coefficient H m becomes sensi- 
tive to the ratio Q/m through its anomalous dimensions which depends on a logarithm of 
v + -n = V- -n = Q/m. Including the RG-summation of these logarithms gives the coefficient 
H m (m, Q I m, fi m , /i) = H m (m, fi m )UH m (Q/m, fi m , fi) appearing in the factorization theorem, 
where Un m is the bHQET current evolution factor discussed below in sections [11 (J| and |TJ. 
Note that in principle H m (m, /x) and the factors in Eq. (E¥) can also have Q/m dependence 
at NNLL. For related discussions see Refs. [|56|, 



Alternatively, the matching coefficient of SCET and HQET jet functions given by H m in 
Eq. ( p4[) can be determined from currents, 

i 1 2 

H m (m, ji) = \C m (m, /i) , (25) 
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where the boosted HQET current is 



Ji(») = C m(m,fi)j£ mET (ii) 



(26) 



with 



bHQET 



(h v+ W n )Y^Y n (Wlh v 



(27) 



The soft Wilson lines Y in this current are the same as those used in the SCET soft function. 
The only distinction is that soft gluons in bHQET no longer couple to massive top-bubbles. 

Due to the large width of the top quarks the B± jet functions can be computed in 
perturbation theory. At tree level they are Breit-Wigner distributions 



B t ™(s,T t ) = Im[B t r(s,r t ) \ =Im 



irm s + iT t 



nm s 2 + r? ' 



(28) 



where we have adopted a normalization such that 



ds Bl ee ' 



The Wilson coefficients in the factorization theorem in Eq. 
at tree level, Hq = 1 and H m = 1. 



(29) 

are also normalized to unity 



B. Factorization of Lifetime Effects and Soft Function Models 



The leading order bHQET Lagrangian is 

C± = h v± [iv± ■ D± — 5m + -T t )h v± . (30) 

In light-cone coordinates, (+, — , _L), we have t>+ = (m/Q, Q/m, 0) and vt = (Q/m, m/Q, 0) 
and gluons/residual momenta scaling as ~ T(m/Q,Q/m, 1) and ~ T(Q/m,m/Q, 1). 
Unlike standard HQET, the ultracollinear gluon fields in bHQET are defined with zero-bin 
subtractions |55| for the soft region. In Eq. ([30]) T t is a Wilson coefficient obtained by 
matching to the full theory and is equal to the top quark total width. This is true to leading 
order in electroweak interactions, to 0(m 2 /Q 2 ) and 0(T/m) in the power counting, and to 
all orders in a s . 3 Finally, 

8m = m po ie — m (31) 



3 Concerning the m/Q expansion this is true because for Q m the hemisphere mass definition is inclusive 
in the top and antitop decay products up to 0(m 2 /Q 2 ) corrections Q. Concerning the T/m expansion 
this is related to the fact that finite lifetime corrections are related to off-shell corrections that are s/m- 
suppressed |5q ]. Concerning the a s expansion this can be seen by carrying out the matching with free 
quark states and noting that the full theory computation of t — > bW gives the total rate. Now only the 
operator of interest (iT t /2)h v h v allows for decays in the effective theory, but it corresponds to a conserved 
current and so does not get renormalized ||. 
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is the residual mass term that fixes the top quark mass definition m that is used in the 
HQET computations. It needs to be consistent with the bHQET power counting [0, 



5m ~ st ~ si ~ T , (32) 

can be computed perturbatively, and is UV- and IR-finite. Note that the way in which 
Eq. (|3"0D will be used is to compute a jet-function where the width smears over a set of states 
of invariant mass mT t ^> Aq CD . Thus, for our analysis there are no AQco/r* corrections to 
Eq. fl30D , just corrections of O(AqcdA^)- 

In Eq. (|2l|) the jet functions B± are expressed in terms of the imaginary part of vacuum 
matrix elements B± in Eq. (0). From C± it is straightforward to see that B± can be 
obtained from the imaginary part of the vacuum matrix element £>± =0 for (fictitious) stable 
top quarks by shifting the energy variable s — > s + % T t , 



B±(s, T u fj) = Im [B±(s, T t , fi) 

= lm[Bi= (s+tT t ,fi)} , (33) 

Here we defined results for stable top quarks, namely the jet function _B± =0 (s,/i) = 
B±(s, 0, /i), and a vacuum matrix element B^r°(s,fi) = B±(s,0, fi). They are related by 

Bi=°(s,n) = Im[Bi=°(s,ti)}, (34) 



and we will refer to i?^ =0 as the stable jet function in what follows. The result in Eq. (|33D 
is in complete analogy to the relation between the production rate of top quark pairs in the 
nonrelativistic threshold region, -E c . m . ~ 2m, where the leading order finite lifetime effects 
can be implemented by the shift -£7 c .m. - > -Bc.m. + i prior to taking the imaginary part of 



the e + e~ — > e + e~ forward scattering matrix element |T6 . 

To separate the different physical effects in the cross section it is convenient to derive a 
factorization theorem for the leading order finite lifetime effects to all orders in a s . To do 
so we define the function 

g{x) = -*-B r ± =°(x,ii) = ~B ± (x,0,m). (35) 

It is analytic everywhere in the complex x-plane, except along the positive real axis, x > 0, 
where the vacuum matrix elements £>^ =0 , defined using Eq. (p2"D with Ft = 0, has a cut 
for intermediate states having invariant masses larger than the top quark mass. Using the 
residue theorem for a contour that envelops the cut, it is then straightforward to derive the 
dispersion relation 

1 r Disc[g(s)] 

9i ~ a) = 2^L dx —^r^ (36) 
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where a is any point in the complex plane not on the positive real axis. With the choice 
a = s + i T t and a change of variable x = s — s', this dispersion relation can be brought into 
the form 

B ± (s,T t ,fi) = [*<& B£*(i-$>,n) I* (37) 

J-oo 71 \ S I 1 t ) 

Note that the upper limit s of the integration can be replaced by +oo since the stable jet 
function only has support for positive values of its energy variable. Equation (|37D states that 
the bHQET jet functions for the physical unstable top quark can be written as a convolution 
of the stable jet functions with a Breit-Wigner function of the width T t . Thus the leading 
order finite lifetime effects can be factorized from the jet function. 4 This means in particular 
that the renormalization properties of the jet functions for stable and unstable top quarks 
are equivalent - a fact that might not be obvious since the evolution of the jet functions 
involves convolutions with distributions. 

Eq. fl37|) reflects the fact that the top quark width acts as an infrared cutoff for the jet 
function through smearing over a Breit-Wigner function In the factorization theorem 



in Eq. (H) additional smearing is provided by the convolution with the soft function, where 
the width of the distribution S(£ + ,£~) is of order the hadronic scale A. Eq. ( pTf ) allows 
us to group both types of smearing into a common infrared function R, with the following 
modified version of the factorization theorem, 

d 2 a _ TT N TT (_ Q 

JmmmI 



(t H Q (Q,fi m ) H m (mj, ,fJm,ti) (38) 

[°°d£ + d£- B r + =° (s t -—, li) B r _=° (s t - 9L } ^ up r r At) . 



The result involves only the stable jet functions and an infrared function defined as 



QkJ Jo J [(^+_/+)2+(2^ t )2] [(£-_£-)2 + (^)2] " 

(39) 

In the r t — > limit we have R(£ + ,£~ ,T t = 0, fi) — S(£ + ,£~,fi), since the Breit-Wigner 
factors reduce to delta-functions. If s^j ^> T t , as in the tail region, then R can be simplified 
with an operator product expansion whose first term depends on the partonic soft function 
S'part which can be computed in perturbation theory. In this region £ ± ~ s t j/Q 3> A, and 
for these momenta 5'(£ ± , / u) = S^tCc , £t) up to power corrections of 0(AQ/s t j). For the 
case T t ^> ^A, R can be computed using an operator product expansion even in the peak 
region, taking Eq. fl39|) with < 1, and again the leading term is determined by .Spaxt. 



4 Note that subleading finite lifetime effects, which are suppressed by r t /m, cannot be factorized as in 
Eq. ( |37| ) since they are not described by a simply shift of the energy into the complex plane []l5|, [l?], |5£| . 
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This is similar to B — > in the multi-scale OPE ]47] , plf where smearing over the soft 
function makes it computable in an OPE. On the other hand, for top-quarks Y t < ^A, 
and the infrared function is significantly effected by nonperturbative effects in 5*. Thus the 
cross-section in the peak region cannot be determined entirely from perturbation theory. 

These properties serve as an important guideline for the construction of a consistent 
model for the soft function to be used beyond the tree-level approximation. They require 
that the perturbative corrections contained in the partonic soft function must be included 
in a viable model in order to obtain the correct leading order term in the operator product 
expansion for the cases £ ± ^> A and Y ^> ^A mentioned above. As discussed in Ref. |3(| 



one way to give a consistent implementation of the partonic soft function in 5* is to use a 
convolution form 

/+oo _ r+oo 
dt / dt S paxt (£ + -£ + ,r-r,v,6) S mod (I + ,t) , (40) 
-oo J — oo 

where 5* part is the soft function computed in perturbation theory at fi, and S mod is a hadronic 
function satisfying 

+oo ^+oo 

d£+l rr d (f + ,t) = i, 

mod /" »+ q — \ ( a \n+m 



oo J — oo 
+oo r+oo 



dt / dt (£ + ) n (r) m S mod (£ + , £') ~ (A QCD )" +m , (41) 

) J — oo 

for n+m > 1. An analogous formula to Eq. ([H]) was used to incorporate moment constraints 
in the study of the soft function in b — > s£ + £~ in Ref. pl|. This form ensures that 5* reduces 
to Spart for £ ± ^> A, and that for all kinematic regions it has the proper fi dependence for 
the MS-scheme. It also gives the proper result for R in taking the limit T 3> — of Eq. (PU[). 
The model in Eq. ( |4"0"D is specified by parameters in S mo( i which involve the hadronic scale A, 
and also by the choice of the scale fi in S part . The convolution generates logarithms of the 
form ln(£ ± //i) and ln(A//i) to be discussed in the next section. There are also complications 
related to removing a u — 1/2 renormalon (as indicated by the subtraction constant 5 in 
S'part) and introducing a renormalon free gap parameter A in the soft function p0| . We will 
use the prescription in Eq. ([40D for our numerical studies of the factorization theorem in 
Sec. [VIII| . Our choice of 5' mod and a review of how the renormalon subtractions work are 



given in Sec. VII 



C. Summation of Large Logs in SCET and HQET 



In this section we discuss the summation of large logs between Q ^> m ^> Y for the peak 
region, and between Q ^> m ^> s in the tail region. We also discuss the ultra-tail region 
where Q ^> s ~ m ^> m 2 /Q. In both SCET and bHQET we can define unitary evolution 
functions C/j, associated with the renormalization group evolution for hard, jet, and soft 
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functions in Eqs. @ and (|T3|). These Ui factors are indicated by the arrows in Figs. ^]and|3|. 
The figures show that there are two ways of doing the renormalization group evolution. In 
the first one, referred to as "top-down", we run the SCET and bHQET production currents, 
starting with matching at a high scale, and running toward the low scales. In the second 
one, referred to as "bottom-up" , we run the individual jet and soft functions, starting with 
initial conditions at the low scales and running up. 

The UV renormalization of the currents in SCET and bHQET generates Uh q and Un m 
respectively. Since the renormalization of a current does not depend on the choice of states, 
these factors do not carry information about the constraints used to define M t and Mj. 
Instead Uh q and Uu m only affect the overall normalization of the invariant mass distribution. 
On the other hand the jet and soft functions have evolution through Uj, Ub, and Us which 
involve convolutions that change their respective shape. In general the jet and soft functions 
also incorporate the prescriptions to define M t and Mj. 5 Hence, although it is expected on 
general grounds, it is not immediately obvious how the running of these functions becomes 
independent of the prescription used to define these invariant masses. The equivalence of 
the top-down and bottom-up approaches ensures that the changes in shape of the jet and 
soft function cancel out in the convolution for any region in /i where they overlap in Fig. |2|, 
and yield the same result as obtained from Uh q and Un m . In the field theory this result 
is encoded in an SCET consistency condition between Uh q , Uj, and U s , and a bHQET 
consistency condition between Uu m , Ub, and U s , where as we will show to all orders in 
perturbation theory the soft function evolution factorizes as 

U s (£+,r) = U s (£ + )U s (£-). (42) 
The consistency conditions between the scales \i and /io are 

\/UH Q (Q,fi,fio) Uj(s, /i , fi) = —U s fi, /i 



/Q \ tti /ins \ 

U Hm y—, ^ Ho) U B (s, /i , fi) = —U s , /i, no) . (43) 

The derivation of Eq. ( f42"D and fl4"3"|) will be given in section [111 C| below, while field theory 
definitions of the Ui are given in sections |111 A| and [111 B . Note that this discussion implies 



that the SCET factorization theorem in Eq. flT3"| ) can be formulated at any scale /i > m, and 
the final factorization theorem in Eq. (|3|) can be formulated at any scale \x < m without 
affecting the renormalization group evolution. For /i > m we have rif — 6 flavors, while for 
/a < m we have nf — 5 flavors for these evolution factors. 

In the peak region factorization theorem in Eq. (§) the B± and S functions are evaluated 
at a common scale /i. Since they involve logarithms of the form ln[(— s — and ln(A//i), 



5 In the class of observables we consider this is the case for S, while J nifi and B± are inclusive because they 
do not depend on the invariant mass prescription which only affects radiation at large angles. 
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FIG. 2: Matching, running, and matrix elements that determine the functions in the factorization 
theorem in Eq.(|J) for the peak region (when s/m ~ T) and for the tail region (when s/m > V). 
The running in Uh q and Un m is local, while that in Uj n , Uj n , Us, and Ub involves convolutions. 
Here the distribution width is T ~ Tt + QA/m. Cases a), b), and c) show three equivalent ways 
to sum large logs with the renormalization group. The consistency equations discussed in the text 
express the equivalence of running from the top-down in case a) and from the bottom-up in case 
b). Case c) is used for our numerical analysis. 

where T t and the hadronic scale A differ, it is natural to consider using different low energy 
scales fir and fi A for the jet and soft functions respectively. In this case we have a region 
between fir and fi A where the consistency conditions no longer apply. Thus for the general 
situation shown in Fig. ||] the factorization theorem becomes 

d 2 a „ ,„ v „ / Q 



dWdM? = a ° Hq{Q ' Hm ( m ' m' /im ' 



(44) 



x / dtd£~ \ d£' + d£'~ Us(£' + ,£'-,fir,fiA) 

J — oo J — oo 

x B + (s t -—,V,fir) B(s- t -—J,fi V ) S(£ + -£' + , r-£'-, fi A ) 
V m / V m / 



or equivalently 
d 2 a 



nr i Mf2 = °o Hq(Q, fi m ) H m ( m, — , fi m} fi A 



Q 



(45) 



x 



x 



/oo 
ds' t ds'i U B+ (s t - s' t , fiA, fir) U B _ (% - s' h fi A , fir) 
-oo 

d£ + di~ B + (s' t - T, fi v ) B„ (s'- t - r, fi T ) S(£ + , £~,fi A ) . 



We will always take fir > A*a (although technically these equations are still valid for the case 
A*a > A*r)- The evolution kernels Ub and Us sum the large logs between fir and fi A , while 
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the large logs that only affect the overall normalization are summed into Hq and H m . In 
Fig. |2| we display three equivalent ways to sum the large logs, labeled cases a), b), and c). 
In case a) we run all terms, from the top-down, from \xq down to /ir, and we run the soft 
function from the bottom- up starting at fi\ and ending at /ir- In case b) we run the soft 
function from fi\ all the way to fiQ, and the jet functions from /ir to fi m , and then from fi m 
to fiQ. Applying the consistency equations between /iQ-/i m and /i m ~/-tr, cases a) and b) are 
equivalent, and both give the result shown in Eq. (JPj). If we take case a) and apply the 
consistency equation between /ir and /iA we obtain another equivalent result, case c), with 
the result shown in Eq. (f4~5|). We will use case c) for our analysis. 

In the previous section we derived a form of the factorization formula fl38f ), which combines 
the finite lifetime effects and the nonperturbative effects into an infrared function R. This 
form gives useful insights for the proper choice of the scales /ir and fi\. In terms of stable 
jet functions and R in Eq. ( |39D the resummed factorization theorem in Eq. (|45|) becomes 

d 2 a ( Q \ 

dM 2 dM 2 = a ° Hq ^ Q) ^ m > Hm l m ' m' /im; ^ A ) ^ 

t t 



x 



/oo 
ds' t ds't U B+ (s t - s t , /i A , /ir) U B _ {& - s' b /i A , /ir) 
-oo 



-oo 

x [Vr Bl=° (s f t -^-, /i r ) B r _=° (s'- t - 9L t f , T t , /i A ) . 

J-oo V mj J \ rrij J 

From the convolution in this result we see that the smearing with R provides important 
information on the infrared cutoff for the fluctuations described by jet functions, and hence 
the choice of /ir that minimizes large logs. Likewise, we see from the definition of R in 
Eq. ( P§D and the form of the soft function in Eq. (HO) that /iA is affected by a smearing 
caused by nonperturbative effects as well as by the scale mT /Q in the Breit-Wigner functions. 
Hence in the peak region we should run down to the scales 

^ (r 1 + ^ + £i£), , A .o( A + ^£ + £M). (47) 

V m m J \ Q Q J 

In principle /ir can be substantially larger than T t depending on the Q/m we are interested 
in. Also with a very large width (which does not apply for the top quark), /iA could be 
substantially larger than the hadronic scale, which would allow for a perturbative prediction 
of the invariant mass distribution in the peak region. For the realistic case of F t ~ 1.5 GeV 
the scale where the logs would be strictly minimized is in the nonperturbative regime, and 
we will specify the soft function at scales fi\ ~ 1 GeV to be close to this regime. In the tail 
region we have s ^ T or s > T and the convolution in Eq. (|H) sets £ ± ~ ms tt t/Q 3> A, so 
to sum the large logs in this region we should instead run down to the scales 



The results above are designed to study situations where s <C m 2 , which is important for 
a precision extraction of the top-mass. In our formalism it is also possible to study the cross 
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FIG. 3: Matching, running, and matrix elements that determine the functions in the factorization 
theorem for the ultra-tail region, in Eq.(|l3|). The running in Uh q is local, while that in Uj n , 
Uj n , and Us involves convolutions. Cases a), b), and c) show three equivalent ways to sum large 
logs with the renormalization group. The consistency equation discussed in the text express the 
equivalence of running from the top-down in case a) and from the bottom- up in case b). 

section in the ultra-tail region, \M t i — mj\ ~ mj, with renormalization group improvement. 
This is the closest analog to the resummation for massless event shapes in regions where 
the jet invariant mass M 2 > QA. In this case we use the SCET factorization theorem in 



Eq. flHf ) which is valid as long as Q 2 3> m 2 , s. Again we take different renormalization scales 
for the jet functions (fi m ) and soft function (//a) as shown in Fig. |||. The SCET factorization 
theorem in the jet mass scheme is [here mj = mj(fi m )] 

d 2 a 



dM 2 dM? 



cr HQ(Q,n m )M(mj } fi r , 



X Jn\St ~ 



• m Ji Hm)Jn{St ~ 



dtdr de ,+ d£'~ u s (e ,+ ,e'-fim,HA) 

J — oo 



(49) 



ds' t ds't Uj{ )Uj{ 

-oo 



X 



~di J„(s' t - Q£ + \mj^ m )Jn{s'i- Qi ,mj, /x m ) S pavt (£ + ,£ ,M,mj) 



With analogy to the discussion above we show three equivalent ways to sum the large logs 
in Fig. H cases a), b), and c). Here cases a) and b) give the first result displayed in Eq. (f49"D 
and are related by applying the consistency equation between \iq and /i m . Case c) gives 
the second result in Eq. fl49|) and is related to case a) by using the consistency condition 
between \x m and /ia- In the region below m where top-bubbles have been integrated out we 
have Uf = 5 active flavors, so rif = 5 for Us in case a), and rif = 5 for U^(fi m , //a) and 
Uj in case c). For HQ(Q,fj, m ) we always have n/ = 6. In the ultra-tail region we generically 
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have s t f ~ "^i 2 , and hence the logs in the jet functions are minimized with /i m ~ m. The 
convolution with the soft function involves momenta £ ~ s/Q ~ m 2 /Q, and hence the 
large logs are summed for /^a — m 2 /Q. 

III. RENORMALIZATION AND ANOMALOUS DIMENSIONS 

In this section we setup our notation and conventions for renormalization of the quantities 
defined in Sec. [II A| in the MS scheme. In QCD the vector current is conserved (although 



in MS one must be careful with the definition f63|), but in the effective theory the currents 
are renormalized. In the running pictured in Fig. |2| we have both local running (anomalous 
dimensions that depend only on conserved kinematic variables) and convolution running 
(anomalous dimensions that depend on variables that can be changed by dynamics in an- 
other sector). Convolution running involves an integration over anomalous dimensions that 
are functions. The coefficients Hq and H m have local running, while the functions J nt n, B±, 
and 5* have convolution running. In this context an example of local running are the loga- 
rithms summed up by the RG-evolution of gauge couplings, and an example of convolution 
running are logs summed by the Altarelli-Parisi evolution equations (which are collinear UV 
logs in SCET). Another important attribute of these functions is whether their anomalous 
dimensions involve ln(/i) factors and hence sum double Sudakov logarithms. These ln(/x) 



factors are induced by cusp angles involving light-like Wilson lines |64|, |65], p6| , 

Thus, in considering the renormalization group evolution in SCET and bHQET the phys- 
ical meaning of the logarithms that are being summed depends on which of four cases we 
are in: 1) local single logs, 2) local double logs, 3) convolution with single logs, 4) convolu- 
tion with double logs. In case 1) we have local running without a ln(fi) in the anomalous 
dimension, and the evolution just corresponds to the change of a coupling constant c(/i) 
from integrating out virtual effects. This is the standard case, well known from the running 
of the gauge couplings and of the electroweak effective Hamiltonian of four-quark operators. 
On the other hand, the UV renormalization in cases 2), 3), 4) are induced by particular 
types of phase space restrictions on real radiation that are built into the effective theory. 
In case 2) we have local running with a ln(/z) in the anomalous dimension, while in case 3) 
and 4) we have convolution running without and with a ln(/z) in the anomalous dimension, 
respectively. These cases are discussed further in Appendix [D|, and will be mentioned as 
they arise in the analysis below. 

A. SCET renormalization 

Top-down running. In SCET we can renormalize the current Jf by switching from a 
bare to renormalized Wilson coefficient, 

C bare = Z c C = C + {Z c - 1)C , (50) 
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where insertions of (Z c — 1)C are treated as counterterms that render insertions of the 
current together with the bare Wilson coefficient UV-finite. Field, coupling, and mass 
renormalization are given by 

= Z l J\ n , A b n ^ = zfA n , m bare = m + 5 m , g h ^ = Z af fg , (51) 

and are all identical to those in QCD H |7|, [70]]. 6 For later convenience we also write the 
mass counterterm as 

5 m = (5 m r le + 5m, (52) 

where (£ m ) pole is the counterterm in the pole-scheme, and for mass schemes other than 
the pole scheme the remainder, 5m = m polc — m, contains finite perturbative corrections. 
Eqs. ( |50"D and fl5"T|) suffice to cancel all UV divergences involving J?. The SCET factoriza- 
tion theorem in Eq. (|H|) is generated by a time-ordered product of two Jf currents. The 
factorization theorem shown in Eq. ( |I3|) only involves renormalized objects. These depend 
on the choice of renormalization scheme in SCET, but this dependence cancels out between 
Hq, J n , J n , and S. The renormalization group equation for C and Hq are 

^■^C(Q,fi) = ^ c (Q,fi)C(Q,n), fi-^-H Q (Q,fi) = 'y HQ (Q,fi)HQ(Q,fi), (53) 

where from Eq. (|50|) 7 C = —Z^ 1 fid/dfi Z c , and since Hq = \C\ 2 we have jh q = 7c + 7c- 
Since the current involves light-like Wilson lines in the n and n direction, the anomalous 
dimension has a ln(/z/Q) cusp anomalous dimension term. The general form is 

1h q (Q, p) = ? Hq [a J In (JL^j + 1Hq [ a J . (54) 

The running of Hq is in case 2) and sums local double logs. Here the current is affected 
by small invariant mass phase space restrictions imposed on real radiation, which leads to 
an incomplete cancellation of real and virtual contributions from soft and collinear effects. 
Once we integrate out virtual effects in the EFT and evolve the current down to the scale 
of these restrictions the cancellation again becomes effective. This is manifest through the 
elimination of large logarithms in EFT matrix elements at the low scale. The process of 
integrating out virtual effects and performing the RG-evolution sums double logs between 
the production scale Q and scale of the phase space restrictions. For the solution to the 
RGE equation for Hq we write 

Hq(Q,h) = H Q (Q,fi Q ) U H JQ,fiQ,fJ.) , (55) 



6 This is true to all orders in a s because there are no zero-bin subtractions ]5q| for the collinear two-point 
functions. To see this note that all soft loop corrections to these functions vanish in Feynman gauge since 
n 2 = 0. Thus there is no region that is double counted and would require a subtraction. 
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where Hq(Q, /iq) is the matching condition at the hard scale of order Q and Uh q (ij>q, /•*) the 
evolution factor with fig > fi. The evolution contained in Ujg is illustrated in Fig. 0. 

Bottom- Up Running. It is well known that there is an alternative but equivalent way to 
renormalize composite operators like Jf , which is often referred to as operator renormal- 
ization (see Ref. J7I] for a review). In this approach the UV-divergences in matrix element 
insertions of the bare operators are absorbed into the renormalization Z-factors multiplying 
UV-finite renormalized operators, (Jf ) bare = ZjJ^. The equivalence of the two approaches 
implies that Zj = (Z~ l ) T ', where the transpose is only relevant in case of a multidimensional 
operator basis. 

Here we consider a variant of operator renormalization that introduces Z-factors for the 
objects J n , J n , and S in the SCET factorization theorem, Eq. (0). We will refer to this 
procedure as factorized operator renormalization. In Sec. [II A| these objects were defined by 
matrix elements of time-ordered products of fields, where each involves a subset of the fields 
contained in the current Jf . To switch from bare to renormalized matrix elements we write 

Jn% e (s) = Jds' Z Jn ^s-s', fi) J n ,n(s',») , 

S bam (£ + ,£~) = Lt'+dt'- Z s (£ + -£' + ,r-£'-,fi) S(i' + ,£'-,fi). (56) 



These equations can be inverted using f ds Zj^(s" — s)Z Jn {s — s') = 5(s" — s') etc. Note 
that the Z-factors only depend on differences of momenta because the renormalization is 
local for position space fields (as discussed further in Appendix |D|). The RGE's read 

d f 

^-^Jn,n{ S ^) = J ds ' ' lJ n ,u{ S ~ S '^) ^n,n(s',/i), (57) 

^S(£ + , f,fi)= fd£' + d£'- 7s(£+-£ /+ , r -£'-, fi)S(£' + , £'~,fi) , 



with the anomalous dimensions being defined as 

7^-^/i) = - Jds" Zjl n (s-s"^)^Zj ntn (s''-s',f,) , (58) 

~fs(t-£' + ,r-£'-,ii) = - fd£" + d£"- Zs\£ + -£" + ,£'-£"-, fi)fi^-Z s (£" + -£' + , £""-£'-, ft) . 

Renormalizability of the theory requires that they are finite as e — > 0, and the general form 
for these anomalous dimension is discussed in Appendix [TJ. For the solutions of the RGE's 
we write 

•/»,»(«, AO = Jds' Uj(s-s',/i,ii ) J n , n {s\ /i ) , (59) 
S(e + ,£-,ii) = fd£' + d£'- U s {£ + -£' + , £--£'-, //, /x ) S(£' + , £'~, /i ) . 
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Note that depending on the set up of scales, as illustrated in Fig. |^, the evolution kernels 
Uj and Us evolve to higher scales (/i > /x ) or to lower scales (// < fi ). 

At any order in a s the anomalous dimensions in Eq. Q59) have the general form 



. 2rfa s 
7 F (t-t = 



t-t' 



+ 1 [a s ]6(t-t') } (60) 



where j is the dimension of the convolution variable t'. Although the soft function anomalous 
dimension has two variables, we will show below in section III CI that it can be written as 



ls {l+, n = S(£ + ) 7 ,(r) + 5{f ) ls (£+) , (61) 



where 7 S has the form in Eq. (p0|). Eq. (pOl) involves a plus-function, which we define by the 
limit in Eq. (|Cl"l) . This is similar to the Altarelli-Parisi kernel in deep-inelastic scattering, 
except for the presence of fd. This explicit dependence on \i must appear to make the plus- 
function dimensionless, and using Eq. (|C3|) can be written as a a ln(/x) factor multiplying a 
bit — t'). Thus they sum double logs, making them fall in case 4), which is a combination of 
case 2) described above and case 3). For case 3) the real and virtual effects cancel for the 
soft contributions, but not for collinear ones, leaving single logarithms to be summed by the 
RGE's. The convolution with plus-functions arises because there are angular restrictions on 
the radiation such as those that occur when an energetic proton state absorbs partons in 
DIS. Summing logs in this case involves a convolution since the logarithms generated by the 
collinear effects depend on a momentum fraction. 

Viewed from the bottom-up each of the jet and soft functions in our factorization the- 
orem has an evolution equation corresponding to case 4). However, in SCET when the 
soft function and collinear-jets are combined in the factorization theorem their convoluted 
product no longer has angular restrictions. So the evolution of the product does not involve 
a convolution. The product still restricts the radiation to small invariant mass and so falls 
into case 2), of local running with double logs as we mentioned above. 



B. bHQET renormalization 



Top-Down Running. Next we discuss the renormalization in bHQET. The renormaliza- 
tion constant for the bHQET current for the counterterm method is defined as 

Cm re = Zc m C m = C m + (Z Cm - l)C m . (62) 

While gluon field and coupling renormalization in HQET and QCD are equivalent, the top 
quark field renormalization differs, with /z.^ are = Z^ 2 h v . The bHQET factorization theorem 
in Eq. (||) is generated by a time-ordered product of two J^hqet currents @ . The soft graphs 
in bHQET are identical to those in SCET up to top-quark vacuum polarization graphs ||, 
and the infrared divergences of the collinear graphs in SCET exactly match those in bHQET. 
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The mass-mode function Ai is IR finite and just enters in the H m matching coefficient. Thus, 
the same cancellation between collinear and soft graphs that yielded local running in SCET 
also occurs in bHQET. So the running of C m is also local. We will demonstrate this explicitly 
in the one-loop computations shown below. 

Next recall that the + and — bHQET sectors are decoupled, so the anomalous dimension 
for C m can only depend on the quantities n ■ v _ = Q/m, n ■ v + = Q/m, and n • n — 2. 
With this theory we are interested in studying small invariant mass fluctuations around the 
top quark mass m. Thus the renormalization group evolution is not related to stronger 
kinematic restrictions on the magnitude of the overall invariant mass of top plus lighter 
degrees of freedom. Here the evolution falls into case 1) rather than case 2). However, the 
anomalous dimension of the bHQET current J^uq ET still contains a remnant of the ln(fi/Q) 
term in Eq. ([54]) in the form of a dependence on \n(m/Q). This /i-independent logarithmic 
term is related to a cusp between Wilson lines. This can be made explicit through the field 
redefinition h v± — > W v± h!£±, where W v± are Wilson lines defined in analogy to Eq. (|23|) and 
hv± are heavy quark fields that no longer couple to gluon fields at leading power. For the 
operator h v+ W n (0) that appears for example in the bHQET current of Eq. ( p?|) this leads 
to J$lw$ + W n (0). Insert ions of this operator lead to the anomalous dimension depending 
on logarithms of the cusp angle n ■ v + = Q/m. [5(| [72|, [73]. Unlike SCET, this angle is 
fixed and independent of \x because the overall invariant mass is ~ m 2 and does not become 
parametrically smaller from the RG evolution. 

The RG equations for C m and H m = \C m \ 2 are 

[i-rC m (m, — , fi) = 7c m (—, fJt) C m (m, — , n 
a// \ m J \m J V m 

fi—H m (m,—,fi) =^ Hm (—,lj)H m (m,— J ij), (63) 

where 7c m = — Z^ 1 fid/d/j, Zc m and 7n m = 7c m +7c m > anc ^ ^ ne general form of the anomalous 
dimension is 

lH m {Q/m, n) = Y Hm [a J In + j Hm [a J . (64) 

We write the solution to Eq. ( |63|) as 

H m (m, — , fj, m , fj) = H m (m, (j, m ) U Hm (—, fhn, lA , (65) 
V m / Vm / 

where H m (m, fj, m ) is the matching condition of the bHQET current at the SCET-bHQET 
matching scale \i m ~ m and UH m {Q/m, /i m , fi) is the evolution factor describing the running 
to a scale fi < \i m - The local evolution generated by Un m is illustrated in Fig. ||. Note 
that the RHS of Eq. (|65|) is not /i m independent at the order one is working, since part of 
this dependence is canceled by the Uq([Aq, fi m ) in Hq(Q, /x m ). This is indicated by the /z m 
argument on the LHS of Eq. (]6"5|) . 
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Bottom- Up Running. Next consider the equivalent approach of factorized operator renor- 
malization in bHQET. In this case we introduce Z-factors for the jet functions B± and the 
soft function S rather than counterterm contributions for the bHQET current. The resulting 
evolution equations for the soft function S agree with those in SCET except for the change 
from rif = 6 to nf = 5, and will not be repeated. To switch from bare to renormalized 
HQET jet functions we write 

B h ± ^(s)= Ids' Z B± {s-s',ii)B±{s',ii), (66) 
where J ds ZZ(s" — s, fi)Z B± (s — §', /i) = 5(s" — s'). The RG equations are 

H-j-B±(s,/j) = Jds' 7 S± (s-s',//) B±(s',fi), (67) 
with anomalous dimension 

/' d 
ds" Z B \ (s-s",fi)fi— Z B± (s" -s',fi). (68) 

The general form for this anomalous dimension can be found in Appendix [D| For the solu- 
tions to the RGE we write 

£±(s,M) = Jd$U B {s-s',n,n T )B±{g,n T ). (69) 

The evolution kernels U B take us from the low-scale fir to a scale /i. Depending on the set 
up of scales, as shown in Fig. Q, we can have \i > fir or /i < ^r- 



C. Consistency Conditions in SCET and bHQET 

In this section we derive the factorization of the soft-function evolution factor in Eq. (f42"D 
and the SCET and bHQET consistency equations quoted above in Eq. ((43|). 

Using Eq. (|56| ) we can obtain a finite result for the SCET factorization theorem by 
determining the UV-divergences for the Z-factors Zj n - and Z$ from each individual SCET 
Feynman diagram contributing to J n S and S. If we instead use the counterterm method 
with the current renormalization factor Z c then a consistent form for the counterterm is 
only obtained once all collinear and soft vertex graphs that contribute to the factorization 
theorem at some order in ot s are added up. Since the two methods render UV-finite results 
and lead to the same predictions, there is a consistency relation between the renormalization 
constants for the operator and the counterterm renormalization method which is very useful 
for practical computations. To derive it we start with Eq. (0) and switch to J^ are , </^ are , and 
gbare us j n g g rs ^ counterterm renormalization and then factorized operator renormalization. 
Equating the results we find that 

\Z C \ 2 S(s-Q£' + )S(s-Q£'-) = Jd£ + d£~ Z^(s-Q£ + ) Z^{s-Q£~) Z^it -£' + Jr -£'') . 

(70) 
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The consistency condition can also be written in terms of the evolution kernels that solve the 
individual RGE's. To derive this form we consider the factorization theorem Eq. ( p~3| ) at the 
scale fj,Q, and use Hq(Q,/i q ) = Hq(Q, ^)Uh q {^, A*o)- Then we write down the factorization 
theorem Eq. (|T^) again at the scale fi and relate the J njfi and S at the scale fi to those 
evaluated at /i using Eqs. (|5TS|). Equating the two results gives the consistency condition 



U HQ (Q,fi,i2 )5(s-Q£' + )5(s-Q£'-) (71) 
fd£ + d£- Uj n (s-Q£ + ,fi, f i )Uj ri (s-Qr,fi,^)Us(e + ~£ ,+ ,£--£'-,fi,^). 



Next we multiply Eq. (|?J) by U J (Q£' ± , /i , //), shift P -> ^ + f ± , and integrate over 
to turn the products of C/j factors on the RHS into delta functions. Carrying out the £ ± 
integrals then leaves 

Us {^ , ij ,fi,fi °) = ^Uvq^Q^i^Uj^s,^,^!/^,^,^ . (72) 
This implies a separable structure for Us to all orders in perturbation theory, so we write 

U S {£ + , £~, /I, fi m ) = U s (£ + , /I, li m ) U s (£~ , /I, /I m ). (73) 

This result for Us implies 



d 



d 



dji 



fi^U s (£ + ,£-)= fi-rU s (£ + ) U s (r) + U s (£ + ) ii—U s {£-) (74) 



d 



d\i 



d£' + d£'- [ ls (£ + -£' + )5(£--£'-) + 5(£ + -£' + h s (£--£'-)]U s (£' + )U s (£'- 



d\i 



so the soft function anomalous dimension has the general form shown in Eq. (|6T|). Now 
using the fact that U Jn = U Jn = Uj, Eqs. (|72|) and (|7J) give the final result for the SCET 
consistency equation 

\/u Hq {Q,h,h ) Uj{s, ho, /i) = —U s /i, ^ . (75) 

This relation expresses the equivalence of running the factorization theorem between /iq and 
/i m from top-down versus from bottom-up as pictured in Fig. |2|. It also states that when 
the convolution RGE's for each of J n , J n , and S are combined as shown in the factorization 
theorem, the result is local running through Uh q without a convolution. This means in 
particular, that the renormalization group evolution of the soft function does not depend 
on the phase space constraints that are imposed dividing up the soft radiation in M t and 
Mi. This is verified explicitly at 0(a s ) in Sec. [IV D| where we show that the anomalous 



dimension of the soft function is unchanged if invariant mass prescriptions are applied that 
differ from the hemisphere prescription. 

Next lets derive the consistency equation in bHQET. Again the use of the factorized 
operator renormalization using Zb ± and Zs correspond to determining the UV-divergences 
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of the individual Feynman diagrams contributing to B + , B_, and S. If we instead use current 
renormalization via Zc m then the consistent form of the counterterm is only obtained once 
all vertex graphs contributing to the factorization theorem at a certain order in a s are 
added up. In analogy to SCET this leads to consistency conditions for the renormalization 
factors and the solutions of the anomalous dimensions. The derivation goes along the same 
lines as in the SCET case, but starting from Eq. (|3|). The consistency condition for the 
renormalization factors is 



5/+ 



c m \ 5[s-^—)5[s-^—) = \dtdt- Zfis- — ) Zb 



Qt 



m / \ m J ./ T \ m / \ m 

x Z s \£ + -£'+,£- -£'-), (76) 



and for the evolution kernels reads 

Q \ *(* Qt' 



^, ft/1 „Wi-^Wl-^) (77) 
d£ + d£~ U B+ (s-^-, fi, fio) U B _ h-^-, n, no) U s (£ + -£' + , f -£'-,», /i A ) • 



Removing the 5-functions and integrals in an analogous manner to what we did for SCET 
above, we obtain the final bHQET consistency condition 



U Hm [— > V> Vo) u b(s, ho, fj) = —U 8 ^— , fi, /i J ■ (78) 

This result expresses the equivalence of running the factorization theorem between \im and 
Ha using either top-down or bottom- up approach, as illustrated in Fig. It also states 
that when the convolution RGE's for each of B + , B_, and S are combined as shown in 
the factorization theorem that the result is local running for H m through Un m without a 
convolution. 

These consistency conditions are important phenomenologically because they state that 
the RG-evolution from the hard scales down to a common low energy scale for jet and soft 
functions does not affect the shape of the invariant mass distributions. Since we have a 
consistency condition in both SCET and bHQET the mass scale m does not affect this 
protection of the invariant mass shape from large log modification. The smooth transi- 
tion between the SCET and bHQET consistency conditions is related to a correspondence 
between geometry and the dimension of the variables in the factorization theorem, as we 
discuss in Appendix 0. Once the B± jet functions reach the scale Ht where their logs are 
minimized, then further evolution of the soft function to ha generates logs that affect the 
shape of the invariant mass distribution. 
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and Ei for the hard, jet and 
and d8lT). Values for the 



EABLE I: Dimension j and anomalous dimensions E[a s ], Eo, 70 ; 
soft functions in SCEE and bHQEE using the notation in Eqs 
one and two- loop cusp anomalous dimensions [74] are shown. For each case our notation for the 
resummation functions u and K is also given. 



D. NLL Resummation 



To sum large logarithms to NLL we must solve Eq. fl54| ) for Uh q , Eq. (|60| ) for U F , and 



Qi 

Eq. ( j64|) for Uu m ■ As discussed in Appendix |B| the general solutions are 
U F (t - t',fi,fi ) 



K 



/'0 



i K (e< E ) 



Q 

m 



e K -> 



m 

~Q 2 



{t-vy- 



(79) 



Here lu = u(fi,fi ), K = K(p,,p, ) and K 7 = K 7 (fi, /i ) are solutions to the integrals in 
Eq. (p8|) . Note that for u>, K, and we use the notation that the first argument \i is 
always the final scale of the evolution while the second argument /^o is always the initial 
scale. This notation is also chosen for the Up evolution factors of the soft and jet functions, 
but differs from our notation for the evolution factors for the current Wilson coefficients, 
Uh and Un m where the opposite ordering is used. Thus writing 



T[a s 



An 



r + 



Ait 



I 2 



Arc 



7o + 



Atx 



1 2 



7i + 



(80) 



we must determine E , Ei, and 70. The determination of E and 70 from one-loop diagrams 
is discussed in sections [TV] and |V| below, and the results are summarized in Table |. To deter- 
mine Ei we make use of the fact that to all orders in perturbation theory F[a s ] is proportional 
to the cusp-anomalous dimension, T[a s ] oc E cusp [a; s ]. The constant of proportionality is also 
determined by the one-loop computations and is summarized in Table ||[ At two-loop order 
the cusp anomalous dimension has the form 



r cusp H 



Air 



-pCUSp 1 

1 



«s(/i) 

Ait . 



1 1 ■ 



(81) 
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with results for ;Tq USP and r™ sp also shown in Table |. For QCD Ca = 3, Cp = 4/3, and we 
have rij flavors. Solving Eq. ( |D8| ) also requires the two-loop /3-function 



where 



H^-a t {ji) = f3[a s 
an 



(82) 



Defining 



Aj = HCU/3 - 2n//3 . 



ft = 34C^/3 - 10C A n f /3 - 2C F n f . 



a s (n) 



a s (fJ>o) ' 

and substituting Eqs. (|80f) and ([82 ) into the integrals in Eq. 
To 



gives 



£*;(//, /x ) 
K 7 (/i,// ) 



l r /W 4vr 1 



lnr . 



j/5o 
7o 



27rr / (r-l-rlnr) 7o /3 ( 



+ '" ' " lnr+ 

a s (/x) 47rr V r o Po 



ri /3A (l-r+lnr) 



47T 



+ 



(83) 



(84) 



(85) 



Info 



In r >. 



Taken together with Eq. (|79|) and the appropriate values for r , 70, and Ti from Table |, 
this determines the NLL evolution kernels. 



IV. SCET RESULTS 

A. Current Matching and Running in SCET 

To determine the Wilson coefficient C of Eq. fllCp we match renormalized QCD and 
SCET S-matrix elements, which we will simply call amplitudes in the following. The QCD 
vertex graphs are given in Fig. [| where momenta p and p are defined. We use dimensional 
regularization for UV divergences and small offshell momenta to regulate the IR divergences, 
letting p 2 — m 2 = p 2 — m 2 = A 2 7^ 0. Since the SCET current should reproduce the infrared 
physics of the QCD current, we can perform the matching with arbitrary external states 
and for any infrared regulator, as long as the same IR regulators are used in the full and 
effective theories. The values obtained for the Wilson coefficients will be independent of the 
choice of IR regulator. 

Results for the QCD graphs in Fig. (| are summarized in Eq. ( |Al|) of Appendix |A|. The 
result for the amplitude includes the vertex graph, MS wavefunction contributions, and the 
residue term, Vfjj + Tf (Z^ — 1) + Tf (R^ — 1), where the subscript 4a on the V indicates that it 
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FIG. 4: One-loop vertex corrections in QCD. 



is the result for Fig. [|a. The form of the residue term is required to ensure consistency with 
physical S-matrix elements. We work in the limit A 2 <C m 2 <C Q 2 . The QCD amplitude is 
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(86) 



where for simplicity we use the shorthand notation A 2 and Q 2 for A 2 + iO and Q 2 + iO, 
respectively. For the SCET computation we have the graphs in Fig. which are evaluated 
in Eqs. (|A5|JA6|) of Appendix |A] with non-zero A 2 = p 2 — m 2 and A 2 = p 2 — m 2 . The sum of 
collinear and soft vertex graphs, wavefunction contribution, and residue is Vfj a + V|b + + 
rf (Z^ - 1) + r<*(i?g - 1). For A = A > and again taking the limit A 2 < m 2 < Q 2 we 
obtain 
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(87) 



7T 



The remaining divergences in Eq. (|8~7j) are canceled by the current counterterm Zq — 1 giving 



1 - 
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4tt 
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(88) 



The running generated by Z c sums In 2 /i terms and falls in case 2) as defined in section 
The renormalized amplitude in SCET then reads 



{P,P\ZcJ? |0)sCET = 
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(89) 



Subtracting Eq. (|8^) from fl86|) all dependence on the IR scales m and A cancels. This is 
an explicit demonstration that at one-loop massive SCET has the same IR structure as in 
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FIG. 5: One-loop vertex and self-energy corrections in massive SCET. Gluons with a line through 
them are collinear, while those without are soft. The soft gluon wavefunction renormalization 
graph vanishes in Feynman gauge and is not shown. 



QCD. Evaluating the difference at the scale \x 
current Wilson coefficient, 



\±q gives the matching condition of the 
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(90) 



As expected from the limit Q ^> m the matching condition is mass-independent and there 
are no large logarithms for fiQ ~ Q. 

The result in Eq. ( p0|) is independent of the choice of the IR regulator and should therefore 



agree with the matching conditions for the massless quark production current. In Ref. [TJ 



the matching coefficient was computed using on-shell massless quarks, and Eq. fl90|) agrees 
with their result. With the regulator used in Ref. |L8| the SCET diagrams are scaleless and 
vanish in dimensional regularization. To see more explicitly how the massless computation 
gives the same matching coefficient we repeat the previous computation with an offshellness 
p 2 3> m 2 , where Q 2 ^> p 2 = p 2 . For this case the renormalized one loop QCD amplitude 
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and from Eqs. (|A5|) and (|A6|) the renormalized amplitude in SCET has the form 
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To obtain Eq. ( p2|) the current counterterm in Zq from Eq. (|88|) was used. Taking the 
difference of Eqs. (|9TD and (0) gives exactly Eq. (0), as expected. 

The imaginary parts in C(Q,hq) and the Z-factor in Eq.(^) arise from real QCD inter- 
mediate states in the QCD vertex diagram that are not accounted for in the corresponding 
SCET diagrams. These SCET graphs account only for fluctuations associated to sectors for 
the n and n directions, while the QCD diagrams do not have such a restriction. Note that 
the complex Z-factor also means that the anomalous dimension jc is complex. However, 
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only | C | 2 appears in the factorization theorem in Eq. (Q) and so the complex phase cancels. 
This treatment is consistent because in the derivation of the factorization theorem the part 
of the phase space integration encoded in the sum over the n and n directions is carried 
out explicitly prior to the formulation of the jet and soft functions in SCET. The matching 
coefficient appearing in the factorization theorem therefore reads 
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To evolve the Wilson coefficient to lower scales we need to solve the RG equation in 
Eq. fl53"[) . The anomalous dimensions are obtained from Z c in Eq. (|88[) and fid/dfj, a s = 
-2ea s + f3(a 8 ), 
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Comparing this result to Eq. (|54| ) we find r o Q = —8C F and 7 = —12C F for the coefficients 
discussed in section [III D| . Also r^ Q [a s ] = — 2r cusp [a s ] and so rf^ Q = — 2I 1 usp . The solution 
for the evolution factor is 

2 

U BQ {Q,n Q ,n)=e*> (^Y°, (95) 

where ojq = Uq(ijl, hq) and Kq = Kq(ijl,}j,q) are determined at NLL order using Eq. fl85[) for 
u lu" and 11 K" . At LL order the solutions are 
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with r = a s {n) / a s {no). Note that solving the RG-equation directly for C(Q,/i) leads to an 
extra phase factor, 

which does not, however, appear in the physical cross section. It's origin is the same as for 
the phase contained in the current matching condition C(Q,//q). 



B. SCET Jet Functions and their Running 

In this section we compute the SCET jet functions J n and J n , defined in Eq. ( |19|) 



perturbatively to 0(a s ). Due to charge conjugation symmetry, the results for J n and J fl 
are identical, so for simplicity we focus on the former. The purpose of the calculation is 
two- fold. First we determine the renormalization factor Zj n , the anomalous dimension r yj n 
and evolution kernel Uj n for the jet function. Second, the renormalized jet function at the 
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FIG. 6: SCET graphs for the one-loop top quark jet function. Dashed lines are n-collinear quarks 
and springs are n-collinear gluons. 



scale fi m ~ m is needed to determine the matching condition of the bHQET jet function, 



which we work out in Sec. |VB| below. Since both running and matching do not depend on 
infrared effects below m we are free to do the computations for stable top quarks. Thus in 
this section we set the electroweak gauge coupling to zero and neglect finite lifetime effects. 

From Eq. ([19]), the tree- level jet function is simply given by the imaginary part of the 
collinear propagator: 



J n (s,m, T = 0,/i) 
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At one loop, the jet functions are given by the imaginary part of the diagrams shown in 
Fig. |6], and results for the individual graphs are summarized in Appendix [A]. We will consider 
the one-loop jet function with and without expanding in s -C m 2 . 

Prior to taking the imaginary part the tree level graph plus the sum of one-loop graphs 
from Eq. ( |A10|) give 
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where s = s + iO and 8m = m pole — m. Hence in the pole mass scheme 5m = 0. The one-loop 
massive jet function J n also appears in the computation for B — > X c iv in the endpoint 
region studied in Ref. [75]. Identifying the combination n + -pu' in Ref. [75] with our variable 
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s we find agreement with their pole- mass result. Expanding Eq. 
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For later convenience we write s = xkj where x is dimensionless and K\ > is a dummy 
scale with dimensions of mass. Taking the imaginary part of Eq. ( P§| ) using the results in 
Appendix y we find that the bare SCET jet function is 
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(101) 



When we expand Eq. ( |101|) for s <C m 2 all terms have a singular (9(s _1 ) behavior except for 
the last two 9{s) terms which are O(s ) and can be dropped in the peak region. The result 
of this expansion agrees with taking the imaginary part of Eq. ( |100| ) . 

To renormalize J^ are the required jet function Z-i actor defined in Eq. ( |56| ) is 
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which gives the anomalous dimension 
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Note that J^ are (s), Zj n (s — s') as well as 7j„(s — s') are all independent of the choice for K\. 
The Z-factor and the anomalous dimension also do not depend on the mass scheme that is 
being employed (determined by 5m). Comparing this result to Eq. (pH) we find Tq"'™ = 8Cf 
and 7 ™'" = GCp for the coefficients discussed in section [T I [TJ|. Also T 
and so r^™'" = 2r™ hp . These coefficients give us the NLL evolution kernel that evolves the 
jet function from the scale /Iq to ji: 
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where uj\ = u>i(fi, //o) and Ki = /Uo) are determined at NLL order from Eq. 

r = a s (/i)/a s (/io). At LL order they are 
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with 

(105) 
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The resummation induced by U j n falls in case 4) , it sums double logs and involves a convo- 
lution. 



Finally, taking into account the counterterm in Eq. ( |102| ) the renormalized jet function 
for a stable top quark at 0(a s ) reads 



J n (s, m, T = 0, fi) = S(s) — 2m 5m 6'(s) + 
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Here x = s/ k\ and the result is independent of the choice of K\ . The last term in Eq. ( |106| ) 
is regular in the small x limit and can be dropped in the peak region. In the ultra-tail region 
discussed in Appendix [F] this term generates the "non-singular" contribution of the SCET 
function. 

From Eq. ( |106j ) we can see that for the variable range s ~ mT further matching and RG- 
evolution is needed for J n : there is no choice of \i that minimizes all the logarithmic terms. 
The particular terms in which the large logarithms appear are controlled by the choice of K\, 
but no choice of K\ removes them completely. For example, with K\ — m and fi = m we still 
have ln(x) ~ ln(r/m); while for = mT and // = K\ we have ln(m 2 //tf) ~ ln(r/m). This 
motivates the matching onto bHQET and RG-evolution between m and V to be carried out 
in section below. For later convenience we quote the leading result for J n when s m 2 
using the choice K\ = m, 
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(107) 



C. Hemisphere Soft Function and its Running 



In this section we determine the 0(a s ) renormalization group evolution of the hemisphere 
soft function, Sh em i(^ + , t~ , fi) and its renormalized partonic expression from one-loop pertur- 
bation theory, S , part (^ + , £~, yu), which is needed to construct the soft function model defined 
in Eq. (0). This model builds in the fact that the full nonperturbative S has the same 
dependence on \i as Spa^t- 

For the computation we use the squared matrix-element expression in Eq. (JT4|) for a no- 
gluon and a single gluon final state. The corresponding Feynman diagrams are shown in 
Fig. [7|, where the solid lines denote the four Y- Wilson lines. Fig. |7|a,b,g are the virtual graphs 
with \X S ) = |0), while Fig. |7|c,d,e,f are the real emission graphs with \X S ) = \e^). Results 
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FIG. 7: Graphs for the hemisphere soft function at one-loop. In this figure the solid lines denote 
Y- Wilson lines, and the line with ticks is the final state cut. 



for the graphs are summarized in Eq. (|A13 ) of Appendix Together with the tree-level 
matrix element the bare hemisphere soft function reads 
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Note that Sp^t £~) is independent of the dummy mass scale k 2 > introduced here. 
However K2 facilitates the application of the standard distribution relation for dimensionless 
variables 
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The relation leads to the following expression for the bare hemisphere soft function: 
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where Gs{£ + ,£ ) contains the finite terms 
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The renormalization factor for the hemisphere soft function then reads 

C F a s r 5{e + -e ,+ )5{£--e'-) 



Z s {£' + -i 
5{t 



+ 0' 



£'--£-) = 5(£ + -£ ,+ )5(r-£'-) 



71 



;ii2) 



e n 



fi 



£+-£'+ 



lxQ{£--£'-) 



while the renormalized soft function is 
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We caution once more that the soft function is in general dominated by non-perturbative 
effects, so the partonic and perturbative result in Eq. (|113 ) can only be used in the framework 
of the soft function model of Eq. ( fKj) or in situations where an operator product expansion 
can be carried out. We also note that to 0(a s ) the renormalized partonic soft function can 
be factored in the form 
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where the partonic soft function with one kinematic variable is 
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The factored form in Eq. (|114 ) was expected at 0(a s ) because the single gluon in the real 
graphs of Fig. [7] can make either £ + or £~ nonzero but not both. Because there is in general 
more than one real parton in the real graphs at 0{a 2 s ) and beyond, the factored form of 
S'part is not expected to hold in general, nor for the full non-perturbative soft function S. 
From Z s in Eq. ( |112| ) one obtains the anomalous dimension 
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As anticipated from the form of the consistency condition discussed in Sec. |111 A| , it has a 
separable structure in the light cone variables £ + and £~ . This separation for 7s(£ + ,£~) 
holds to all orders in ot s as discussed in section [111 C\ Comparing this result to Eq. (^) we 
find Tq = — 4C F , 7q = 0, and infer = — r™ sp for the coefficients discussed in section [111 D . 
The anomalous dimension of various soft-functions in SCET were studied in Ref. [|7E[], and 
in particular the one-loop anomalous dimension for the jet-energy soft-function in e + e _ — > 
dijets was derived. This anomalous dimension has the same form but opposite sign of the 
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anomalous dimension in Eq. ( |116| ). However this is inconsequential, since there is no simple 
relation between the jet energy soft function and the hemisphere soft function studied here. 
The hemisphere soft function's 7 S is the same as the anomalous dimension for the soft 
function for Drell-Yan in the endpoint region at one-loop H3j. 

Before solving the soft-function anomalous dimension we pause to consider the number 
of active flavors in the RGE. For the SCET computation we can consider the five lightest 
flavors to be massless. In addition because we are above the scale of the top-mass, we have 
a top-quark loop contribution to the soft-gluon induced /3-function. The top-quark bubble 
couples to the soft-gluon with a multipole expansion since the soft-gluon has p 2 <C m 2 
for the purpose of power counting. Thus the bubble enters as an insertion of the vacuum 
polarization function at zero-momentum on the gluon line, 11(0). The renormalization of 
this n(0) means that for /i > m there are a total of rif = 6 flavors for the SCET running 
of the soft-function. For fi < m the contributions from these top-bubbles is integrated out, 
and the soft-function runs with nj = 5 flavors. Thus, in particular we always have rif = 5 
for the soft-function in bHQET. 

To solve the RG-equation in Eq. fl5T|) we can use that the same equation holds for 
the evolution kernel Us(£ + , £~, /i, /io) defined in Eq. (|59|) that describes the running of 



the soft function from /i to the scale /i. Using the separable form Us{£ + 
U s (£ + , /i, fio)U s (£~ , /i, /io) of Eq. (fT3|) one obtains from Eq. ( |116| ) the relations 
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Here K is a separation constant that can be set to zero. 7 The solution for U s (l , /i, /m) is 
given by Eq. (ff|), 
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where at NLL order we use Eq. ( p5|) for u> 2 = u 2 ([i, fio) and K 2 = K 2 (n, /io) with the values 
of r 0i i and 71 determined above. At LL order these are 
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where r = a s (fi) /a s (fi ). The running generated by Eq. (|118|) falls in case 4). Note that for 
the jet function a>i(/i, /io) > for /i > /j , while for the soft function u 2 (fi, /io) < 0. Although 



7 Note that keeping the ±K term simply adds a multiplicative factor of (h//ia) ±K to the solutions, and 
cancels in the product U s (£ + ,£-) = U S (£+)U S (£-). 
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this affects the behavior of the jet and the soft functions for large values of their arguments, 
the convolution of the jet functions and soft function always remains finite. 

The factorization of the soft function evolution, Us(£ + ,£~) = U S (£ + )U S (£~), is necessary 
for the consistency equation in Eq. flTo] ) to hold since it allows the cancellation to indepen- 
dently occur for the two jets, which are each convoluted with one of the variables of the soft 
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This verifies that the SCET consistency condition is satisfied for the NLL evolution factors. 
The consistency equations can also be verified at the level of the distributions in Eq. (ffll) 
using the results in Appendix [D| Between \iq = 5m and \i m = m = 172 GeV we find 
< co>o(mQ)M) ~ 0.14 at LL and NLL order. 



D. Universal Running for a Class of M t f Definitions 

In the last section we showed that the soft function for hemisphere invariant masses 
satisfies the SCET consistency condition. Since the renormalization is not sensitive to low 
energy properties of the soft function, like the mass definitions, one should expect that there 
is a broader class of soft functions that are consistent with the RGE in our factorization 
theorem. In this section we demonstrate this explicitly by working with a broader set of 
mass definitions and calculating the corresponding soft functions to 0(a s ). 

In Ref. it was shown that the form of the factorization theorem in Eqs. ([|) and (p5[) 
is retained for any M t> i prescription that assigns the hard top and antitop decay jets un- 
ambiguously to M t and Mj, and the momentum of every soft particle to either M t or Mj. 
The former condition ensures that the jet functions B±, which are fully inclusive for the top 
decay products and collinear radiation, remain unchanged. The latter condition ensures that 
concentrating on M t) j in the peak region automatically selects events in the dijet region for 
which the SCET-bHQET setup can be applied. The first condition is satisfied by reconstruc- 
tion methods since for Q ^> m the hard jets and collinear radiation are collimated to two 
back-to-back regions of the detector and the decay products only have a power-suppressed 
probability of 0(m 2 /Q 2 ) to show up in the opposite hemisphere 0. The second condition 
restricts us to invariant mass definitions that incorporate all soft radiation. Examples in- 
clude the hemisphere definition used in the last section, and particle recombination methods 
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such as those based on kr jet algorithms JF7J with a y cu t parameter chosen so that all soft 
radiation is assigned to the hard jets from the top/antitop decay. Based on the equivalence 
of the top-down and the bottom-up approach to the renormalization of quantities in the 
factorization theorem, and the fact that the renormalization of the top-antitop production 
currents can only depend on virtual corrections, it was concluded in Ref. that the renor- 
malization properties of this class of soft functions does not depend on the prescription how 
the soft gluon momenta are assigned to M t and Mi. 

Lets now extend the soft-function analysis beyond hemisphere masses, by setting up a 
more general definition for M t and Mj and hence for the matrix element defining S(£ + , E~, /x). 
Since the contributions from the virtual graphs in Fig. |7]a,b,g are unaffected by phase space 
constraints it is sufficient to consider the graphs [?]c,d,e,f describing real soft gluon final 
states. It is useful to write the gluon phase space integral given in Eq. (|A12|) in terms of the 
perp-momentum q± and the angular variable 

Q 

x = tan- = e"\ (122) 

where 9 is the gluon angle and rj the rapidity with respect to the top momentum direction. 
This gives 
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where fl is given in terms of /i in Eq. (K3j). For the hemisphere invariant mass prescription 
gluons in hemisphere-a (0 < x < 1) are assigned to the top, and gluons in hemisphere-b 
(1 < x < oo) are assigned to the antitop. One can interpret this hemisphere prescription 
as a crude jet algorithm. A prescription such as the k T jet algorithm [|77fl , that is tuned 
such that the total number of final jets equals the number of hard jets from the top and 
antitop quark decays, leads to a more complicated pattern since it depends on the particular 
momentum configuration of the hard jets. 

However, the situation is simplified since at leading order in the power counting the 
hard jets are assigned unambiguously to the top and antitop invariant masses. Thus upon 
averaging over all hard jet configurations the jet algorithm assigns a soft gluon to either M t 
or Mi according to a probability function, f(x), that depends only on the angle 9. Expressing 
the phase space integration of the 0(a s ) soft function for this general jet algorithm we have 
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where f(x) gives the probability that a soft gluon with x is assigned to M t . For the hemi- 
sphere masses we have f(x) = 0(1 — x). Consistency at leading order in the m/Q power 
counting requires that /(O) = 1 and /(oo) = 0, i.e. the soft gluon is assigned with unit 
probability to M t (Mj) if it is radiated in exactly the top (antitop) momentum direction. 
Using the identity of Eq. (|109|) and the scaling variable K2 from the previous subsection 



it is then straightforward to determine the bare soft function: 
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where it is the finite terms that depend on the arbitrary probability function f(x) 
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For the hemisphere masses f n = for any n, and Gs(£ + ,£~) reduces to Gs(£ + ,£~) in 
Eq. Qlll]) . For a more general prescription for the soft gluon assignments that is symmetric 
under exchange of top and antitop, i.e. has f(x) = 1 — f(l/x), then one still has fo = 0, 
while f\ is in general non-vanishing. 

The result in Eq. ( |125|) demonstrates that the UV-divergences and the RG-evolution of 
the soft function are not affected by the phase space constraints imposed on the soft gluons, 
whereas the UV-finite contributions depend on them. This demonstrates to 0(a s ) that the 
form of the factorization theorem in Eqs. (§,^) is retained for the class of invariant mass 
definitions described above, and that different mass prescriptions only affect the form of the 
soft function, but not its renormalization scale-dependence. 

V. HQET RESULTS 



To describe scales below the top mass we need to integrate out m by switching from 
SCET to bHQET. Here we describe the bHQET analogs of the matching, running, and 
matrix element results given in the previous section on SCET. 
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FIG. 8: Nonzero one-loop vertex and wavefunction corrections in boosted HQET. Graphs a), b), 
c) involve heavy quark fields h v± , while graphs d) , e) , and f ) only involves the Wilson lines Yn and 



A. bHQET Current Matching and Running 

In this section we determine the matching and the running of the ti current in bHQET 
at 0(a s ). For this we need to consider the one-loop graphs in Fig. |[ For convenience, the 
relevant bHQET Feynman rules have been collected in Appendix ||, as are the results for 
the individual graphs. We use dimensional regularization for UV divergences and offshell 
momenta to regulate the IR divergences. For the top and antitop quark momenta we take 
p = -i- r^. nnrl fjf* = mv^ + rl, respectively, and then let p 2 — m 2 = 2mv+ ■ r+ = A 2 , 



p 

and p 



mv+ + r+ and p M 

>^2 



m 



m 



2mv- 



A 2 , with A ^ 0. 



The sum of the three vertex contributions, the wavefunction contribution and the residue 
is Via + Vjb + Vic + T^(Zh — 1) + rf (Rh — 1) and gives the bare amplitude 



(p,p|J^|°)bHQET 



ln ; 



. a s C F (2 m 2 2 
^ + 4ln 2 ^ + 4ln^ + 44 



;i28) 



where Q 2 = Q 2 + iO and A 2 = A 2 + iO. The UV divergences in the bHQET current are 
subtracted by the counterterm for the Wilson coefficient 
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giving the renormalized bHQET amplitude 
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As discussed already in Sec. [I I [B| the current renormalization constant contains a term 
\n(m/Q)/e with a coefficient that agrees with the coefficient of the \n(fi/Q)/e term in 
the renormalization constant of the SCET current in Eq. fl88|). Also, the anomalous di- 
mension of the bHQET current only exhibits single 1/e poles and thus sums only single 
ln(/i) contributions. Although this running sums double logarithmic terms of the form 
\n(m/Q) ln(yu/m)] fc , it formally belongs to class 1). 
From the difference of the renormalized bHQET amplitude, (p,p\Zc m Jt\fy ■* an d the 
renormalized SCET amplitude in Eq. (|85|), we obtain the bHQET current matching condi- 
tions at the scale /i m : 



C m (m,n r 
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The matching coefficient H m (m,fi r , 
orem reads 
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= \C m (m, fi m )\ 2 that appears in the factorization the- 
a s C F 
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m 2 m 2 



71 
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(132) 



This matching result only depends on the parameter m, and at the scale /i m ~ m there are 
no large logarithms in H m (m, /i m ). 

The anomalous dimension is obtained from Zc m using Eq. (BHT) and gives 



a s C F 
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(133) 



Comparing this result to Eq. (Q) we find r o m = — 8CV, 7(f m = —8Cf, and infer 
_2pcusp £ or ^ e coe fg c j en ts discussed in section |111 D\ The solution for the evolution factor 
for the mass scale coefficient H m in Eq. (^) reads 



U, 



Q 



m 



,K W 



Q 2 



where at NLL order we use the expressions in Eq. 
K 00 (fj,, jj, m ). At LL order we have 



(134) 

for u = uj (fi,fi m ) and K 00 = 



Ad 



In 



q g (/i) 

Oi s (/^rn) 



(135) 



and just as in the running with Uh q , Uq L ([1, Av) = (^Cp/Po) ln[a s (//)/a! s (// m )]. Note that 
as in the case of the SCET current the RGE solution for the current Wilson coefficient 
C m (m, fiQ, fi) contains an extra phase factor, 



C m (m,n) = yjH m {m,n) 



(136) 
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FIG. 9: bHQET graphs for the top quark jet function. 



that does not, however, appear in the cross section. 8 The origin of this phase, and the 
reason it drops out of the final predictions, is the same as for the SCET current Wilson 
coefficient discussed in Sec. [IV A| . 



B. bHQET Jet functions Matching and Running 

In this subsection we determine the bHQET jet functions B± defined in Eq. ( pip at 0(a s ), 
obtaining one- loop corrections to the Breit-Wigner distributions in Eq. (|28|). The results for 
B + and B- are identical by charge conjugation. We also determine the bHQET jet function 
renormalization factor the jet anomalous dimension 75, the NLL evolution kernel Ub, 
and finally B± at NLL order. By comparing the jet functions in bHQET and SCET we 
confirm that their IR divergences agree. Finally, we demonstrate that the matching condition 
for H m (m, fig), already obtained for the top-antitop currents in Eq. ( |132| ), is reproduced by 
jet function matching. This is a reflection of the statement that we have the same soft 
function in the SCET and bHQET theories to 0(a s ). Thus the soft function computations 
in sections [IV C| and |IVD| apply equally well for bHQET. 



For the computation it is convenient to use the formulae from Sec. [II B| which determine 
the jet function for the unstable top quark from the results for a stable bHQET theory. To 
do this one can either use the relation Eq. ([33]) which shifts the invariant mass variable into 
the complex plane, or use the convolution relation in Eq. (|37j). 

The bHQET jet functions are given by the imaginary part of the vacuum matrix elements 
B± defined in Eq. (^). At tree level they are just given by the the HQET propagator, 

B ± (s,V t = 0) = -— — L-. (137) 
7rm s + 2O 

At one loop the diagrams contributing to the vacuum matrix elements B± are shown in 
Fig. (§). Results for individual graphs are given in the appendix. The sum of the one- loop 



It is interesting to note that the result in Eq. (|136| ) can be obtained from running the heavy-to-heavy 
current in HQET S, analytically continuing to the production region |7q], and expanding in m/Q. 
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graphs in an arbitrary mass scheme gives the bare expression 
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(138) 



In general the residual mass term 5m = m po i e — m is nonzero and uniquely fixes the 
mass scheme m that is being employed in HQET. In an arbitrary mass scheme we have 
B±(s,T t , fi,5m) = B±(s — 25m,T t , ft). Here 5m is computed as a perturbative series in a s . 
We have 5m ~ a s at lowest order, and the 5m in the HQET Lagrangian, Eq. (|30|), should 
be included as a perturbative insertion. This yields the term shown in Eq. ( |138j ). 

The result in Eq. ( |138[ ) can be compared to the computation of initial state radiation 
from a heavy color scalar resonance produced by the collision of massless colored and neutral 
scalars in Ref. [HJ. At leading order in the 1/m expansion the HQET gluon interactions 
are spin independent and so only affect the normalization. Furthermore to 0(a s ) there is 
no difference between initial state and final state radiation, the signs of the iO terms in the 
eikonal propagators do not modify the result. In our calculation the analog of the initial 



state colored scalar in Ref. [15| is given by the final state Wilson lines in our jet-function. 
Thus, we expect that the linear combination of terms in Eq. ( T38|) in the pole mass scheme 
where 5m = should be the same as obtained in the scalar computation |fl5[| , and we have 



checked that this is indeed the case. The scalar analysis of Ref. was based on diagrams 
rather than deriving a factorization theorem, so an operator analogous to the one we give 
for our jet-functions was not given. 

The renormalization of the vacuum matrix element and the jet functions for the stable 
or the unstable bHQET theory is equivalent, so one can obtain the renormalization factor 
for the jet functions from Eq. (|138| ). The result reads 



Z B± (s - s') = 6(B - s>) + ^{ 5{s - s') 
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(139) 



Note that the bHQET jet functions Z-factor and their anomalous dimension do not depend 
on the mass scheme that is being used. The renormalized vacuum matrix element has the 
form 



B±(s, 0, fj,, 5m) 
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(140) 



The renormalized jet function accounting for the large top quark width is then either given 
by taking the imaginary part of Eq. ( |140[ ) upon the shift s — > s + iT t (see Eq. (|3~3"D) or by 
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applying the convolutions formula of Eq. (|37| ) . The result is 
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where the arctan is evaluated in the third quadrant. In practice, in the presence of the width 
it is more convenient to not bother evaluating the imaginary part explicitly, and simply use 
B±(s, T t , /i, 6m) = lm[B±(s+i T t , 0, /i, 5m)~\ . We also present the stable bHQET jet function, 



B± °(s, ji, 5m) = Im [B±(s, 0, //, 5m)] 
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where we have allowed for an arbitrary rescaling of s = K3 z. A convenient choice for the 
parameter K3 is K3 = /1, where z = s//i and 
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To determine H m we can match the bHQET and SCET jet function results for s = 
s/m <C m. For the matching the top width only takes the role of an IR parameter, and 
the computation is most conveniently carried out for stable top quarks. If the same mass 
scheme is used in SCET and bHQET, then the 5m terms are the same, and cancel in the 
matching. The jet-function matching coefficient can be obtained from matching either the 
jet functions or the vacuum matrix elements. For this computation it is convenient to pick 
k 3 = m in Eq. ( |142| ) and subtract it from Eq. ( |107| ) to obtain 
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Using H m = T + T- this agrees with Eq. ( |132| ) at 0(a s ). 

The anomalous dimension for the jet function is determined from Eq. ( |139| ) and reads 
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Comparing this result to Eq. (|60| ) we find Tq = ACp, Tjf = ACp, and infer Tf = r™ sp for 
the coefficients discussed in section |111 D| . The solution for the evolution equation ( j69| ) is 



U B (s - s',/i,/i ) 
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where lo\ = ui(fi, /i ) and K s = K 3 (/i, /i ) are determined at NLL order using Eq. 
r = a s (fj,)/a s (fi ). At LL order 
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In comparing the bHQET jet function to that in SCET, the most striking difference is 
that the dimension-1 variable s is natural for bHQET, whereas in SCET we had a natural 
dimension-2 variable s. This causes a difference in the convolution of jet and soft functions 
in the two theories. However, comparing the bHQET jet function evolution function Ub to 
the SCET evolution function Uj in Eq. (|104|) , one notices that we have the same function uj\ 
of a s . This is crucial to the fact that the spectrum remains protected against large logs as we 
evolve below \l = m, with the consistency condition in SCET carrying over to a consistency 
condition in bHQET, as given in Eq. ([Ef). In particular, the relation between the soft and jet 
evolution factors, uji(li , fi) = u; 2 (/z, /z ), remains valid. This ensures that the plus functions 
in Uj and U s match, which was the key to verifying the SCET consistency condition in 
Eq. (|121|) above. The fact that u\ is unchanged can be viewed as a compensation between 
a change in the geometry of the physical color flow governing the QCD dynamics in the 
two theories, and a change in the dimension j of the natural jet function variable used in 
Appendix [TJ. Here the geometric properties are encoded in the cusp anomalous dimensions 
T , which are determined by the cusp angle between Wilson lines in the jet functions. In 
particular the equality of the uj\s follows from the equality of the ratio 
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as can be seen from Eq. fl55| ) in Sec. |111JJ| . To verify the bHQET consistency condition to 
NLL order we use ujo(li , li) = — u;i(/i ,/i) = — ^(a*, A*o) an d note that 
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which allows us to obtain the desired result in Eq. 
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Again the equivalent bHQET consistency equation with distributions given in Eq. (^Ffj) can 
be verified using results from Appendix 0. Between Li m = 172 GeV and li — 1 GeV we find 
< u}q(liq, fx) < 0.93 at LL and NLL order. Interestingly, from the factorization theorem in 
Eqs. fl44D and ( fl"5|) the Ub factor is run down from /ir to fi\ giving lo\ < 0, or the Us factor 
is run up from li a to fir giving u 2 < 0. Thus we never exceed the bound uj < 1 on the range 



of validity of the convolution resummation formulas given in Sec. [II D 
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Finally, we quote analytic results for the resummed bHQET jet functions. Using the tree 
level result for bHQET propagator B± in Eq. ( |137| ) as the initial condition at the scale // it 
is straightforward to determine the LL result at the scale // by carrying out the integral in 
Eq. (|5§D. For the vacuum matrix element this gives 
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where ujY" and K^ L are given in Eq. ( |147| ). Now using Eq. (|37|) gives the LL bHQET jet 
function at the scale /i with the tree level jet function B± ee = T /[nm(s 2 + T 2 )] as the input 
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Since the boundary condition is specified at tree level the result does not involve 5m. When 
the 0(a s ) jet function is taken as the initial condition at the scale /io with NLL evolution 
we obtain what we will call the NLL result for the jet function. This result depends on 5m 
and is given by the analytic result 
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where H(u) is the harmonic number function and \fr'(z) = d/dz[T'(z)/T(z)] is the derivative 
of the polygamma function. This result was derived using Eq. ( |E6|) from Appendix [E|. 



VI. SHORT-DISTANCE TOP JET MASS AND THE TOP QUARK POLE 

One of the main goals of the factorization based analysis of d 2 a/dM 2 dM 2 is to facilitate 
a high precision determination of the top-quark mass. To do so it is important to explore 
the correspondence between the top-mass and the peak in the predicted M^f invariant mass 
distribution. Eq. (H) shows that the peak position of the top/antitop invariant mass dis- 
tributions is affected by perturbative corrections in the jet functions B±(s, T t , //) and by 
nonperturbative effects through the convolution with the soft function S(£ + , £~, //). In this 
section we analyze effects of the jet-functions on the peak, and define a consistent short 
distance jet-mass scheme. The jet functions were computed to 0(a s ) in section [V[ In these 
computations the top quark width Ft provides an IR-cutoff that makes the perturbative 
determination of the jet function line shape valid for any value of s, and in particular for 
the peak region. 
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As can be seen from Fig. ^ the Feynman diagrams for B± contain the HQET heavy 
quark self-energy. The other diagrams are generated from the Wilson lines W n ^ and render 
the jet functions gauge-invariant. In the pole mass scheme the HQET self-energy develops a 
linear sensitivity to IR momenta |7Dj. This is caused by an ambiguity of O(Aqcd) in the pole 



■m 
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c)/ m 



polo 



mass, and hence in the pole-scheme invariant mass variable s, i.e. in s = (M t 2 . 
where 8m = 0. In perturbation theory this ambiguity is associated to an asymptotic behavior 
oc /ia s (/i) n+1 /9o n\, where /x is the scale employed for a s (p) in the jet function. The situation 
is similar to the total cross section for top-antitop pair production in the threshold region 
Q ~ 2m ||], |81| , where using the pole mass the peak position can not be rendered stable 
in perturbation theory. Thus, in the pole mass scheme B± is expected to have a poorly 
behaved perturbation series, indicating that it is not the pole mass that can be accurately 
determined from the measured invariant mass distribution. This is because the pole mass 
is defined order by order to be a zero of the inverse heavy quark two-point function that is 
not observable physically. This feature is demonstrated for the one-loop jet function below 
and will be analyzed at higher orders in Refs. [^, |83|j . 

It is therefore advantageous and even mandatory to switch to a short- distance mass 
scheme that can stabilize the location of the jet function peak location in perturbation 
theory. In a general mass scheme the location of the peak is determined by 



dB±(s, T t , /x, 8m) 
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At tree level 8m = and the jet functions are equal to the Breit-Wigner functions in Eq. (p8|) , 
so Sp Ca k = 0. At 0(a> s ) the jet functions B± are given by the expressions in Eqs. ( |141|) and 



solving the condition in Eq. ( 154j) perturbatively gives 
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As explained in Ref. a viable short- distance mass scheme must have 8m ~ a s T t in order 
not to violate the power counting. This condition rules out the MS-mass as a candidate since 
^ m Hs ~ a s m ti and hence violates the factorization theorem for the cross-section. A viable 
short-distance mass scheme can be defined using Eq. ( |154|) by demanding that s pca k = 
order-by-order in perturbation theory. Eq. ( |155| ) then determines 8m at NLO. This mass 
satisfies the power counting criteria, but unfortunately has a complicated dependence on the 
renormalization scale. This can be seen from the peak position at LL order, derived using 
Eq. ( |152| ) for B± with the initial condition of a Breit-Wigner at the scale /Xq. This LL result 
is independent of the mass scheme and we find 



Vak (^o) = ~ r * COt 



TT 



a s (/x) 



Po La s (/x )j 



+ 



(156) 



L2 + ^ 1 LL (/x,/x ). 

Here (/x, /x ) is given in Eq. ( |1471 ), and in the second equality we show the leading term 
for small uii. Due to the non-linear nature of the cotangent in Eq. fll56| ) a mass scheme 
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determined by s pea k = is not transitive, in the sense that s^ k (fi, /ii) + ^^.(//i, n?) ^ 
Spg ak (/i, yU 2 ). This makes a mass definition based on the peak position awkward to use. The 
problem occurs because the peak position is a local feature of B±(s,[i), while B± requires 
a convolution for its RG-evolution as seen in Eq. (|69| ) . If Eq. ( pL54j ) is evolved to a different 
renormalization scale then it involves an integral over B±. In perturbation theory this 
nonlocal feature is reflected by terms [CVa s ln] fc in the expansion of Sp^ ak . 

To define a transitive short-distance jet-mass mj(/i), which is still closely related to the 
peak-position, we will use the first moment of 5^ =0 . To LL order it suffices to simply use 
an upper cutoff L m on this moment, and define 5m j so that this moment vanishes 
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Different choices of L m define different schemes for the mass. As indicated, it also suffices 
to define the mass scheme using the zero-width jet function. As shown in Eq. fl37|) the jet 
function for a non-zero width is related to the stable one by 
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and so the stability of B± is directly transferred to B±. We can solve Eq. 
only the linear term in 5m j, thus 
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Integrating by parts gives the solution 
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[2 ds Bl=°(s, 0)] - 2L m Bl=°(L m , 0) 
Expanding in a s to one-loop order we find 
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To obtain a consistent mass for top-quark jets we must choose the scheme parameter L m ~ 
T t . We will adopt L m = 1 GeV for our analysis since then 5m j gives |Sp^| < 32MeV 
for = 2-10 GeV and hence a very stable peak position. The use of a first moment to 
define a mass scheme as in Eq. ( |157 ) has been applied in a similar way earlier to inclusive B- 
decays f83|, to give what is known as the shape-function scheme. However the shape- function 
scheme does result in a different mass definition from the jet-scheme defined here. 

To achieve cancellation of the O(Aqcd) renormalon ambiguity in the jet mass scheme the 
scale \i in 5mj(fi) needs to agree with the renormalization scale used for the strong coupling 
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FIG. 10: (a) The bHQET jet function mB±(s,Tt, n) as a function of s at tree level (black long 
dashed line) and 0(a s ) in the pole mass scheme (green short dashed lines) and the jet mass scheme 
(red solid lines) for \x = 2 (lower lines) and 5 GeV (upper lines), (b) Imaginary part of mB±(s, 0, /x) 
in the pole mass scheme for fi = 2 GeV plotted in the complex s-plane. Solid green lines indicate 
Re(s) = or Im(s) = in the plane where Im(mB±) = 0. For the strong coupling we used 
a s = 0.262, 0.203 for (j, = 2, 5 GeV. 



in the corrections of the jet function. Using Eq. ( pi]) the one-loop relation between the pole 
and jet mass rrij is: 
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We can also derive a LL result for the running jet-mass. From the NLO 5mj in Eq. (|161|) 
we can compute a renormalization group equation for mj(fi), whose LL solution is 
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To verify that this result contains all the leading-logs, we use Eq. O to determine (5£ =0 ) LL 
with evolution from up to fi. Using the LL jet function in Eq. ( |152| ) and solving Eq.( |157| ) 
without expanding in 5m, we find a solution Sm(fj,, fii) that contains all leading- logs between 
Hi and ii. Then by taking mJ L (/i) = mj(/i ) + 5m(fj, , fj,i) — 8m(fi, //i) we obtain a 
independent result that reproduces exactly Eq. (|163[) . The jet-mass mj{p) has a standard 
series of [/9oQ!s m ] fc terms, and as far as its RG-evolution is concerned behaves very similar 
to an MS mass. In particular, this jet-mass is transitive at LL order. 

Results for the jet functions are shown in Fig. pTj|a, where we have plotted mB±(s, T t , //) 
at tree-level (long dashed black line), and at NLO in the pole mass scheme (green short 
dashed lines) and in the jet-mass scheme (solid red lines). We take T t = 1.43 GeV. At 
0(a s ) we use Eq. ( |141|) , taking 5m = in the pole mass scheme, and taking 5m j with 
L m = 1 GeV from Eq. ( |161| ) in the jet-mass scheme. For each 0(a s ) prediction we show 
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two curves, one for fi = 2 GeV (lower lines) and one for fi = 5 GeV (upper lines). While the 
resonance peak is located at s = at tree-level, in the pole scheme at one-loop it is shifted 
by 250 MeV towards smaller masses. In the jet-mass scheme the peak is located at s ~ 0. 

One may wonder how the shift of the jet function in the pole scheme arises, given that 
B± in Eq. (141) obviously has a pole at s + iT t = 0. The reason this pole is not visible in 



Fig. |IO| a is that jet-function is modified by powers of ln[(— s — iY t )/ fi\. To illustrate this, 
consider the inverse of the stable vacuum matrix element in the pole mass scheme 



— I , ... L CX s Cf 



[B±(s, 0, fi)] = -ixm (s+iO) < 1 



47T 



41nM-/M+41nf-A-U4+ 57r2 



iOJ \—s — iO/ 6 



(164) 



In Fig. pT0| b the imaginary part of mB± is plotted in the complex s-plane for fi — 2 GeV. 
The small positive peak visible at s = is related to the zero of (£>±) _1 at s = and 
thus connected to the pole mass. However, it is inaccessible physically when the finite top 
quark width is accounted for, i.e. when B± is evaluated in the upper complex half-plane at 
s + iT t . Moreover, it has the wrong causality structure since it leads to a small negative 
dip when approached from the upper complex half-plane. The vacuum matrix element is 
instead dominated by the pole on the negative real axis, that is visible as a large peak to 
the left of the smaller peak at s — 0. Conceptually this means that the physical pole of the 
jet function is not located at the pole mass and that the pole mass per se is not tied to a 
physical object. This conclusion is fully compatible with previous work on the consequences 



of the pole mass renormalon problem fl85|, [Sq, 87, BR 89. 90], from which it is already known 



that the pole mass is an unphysical parameter. The analysis of the calculable B± provides 
a surprisingly direct view on the mechanism of how this is achieved within perturbation 
theory. 

From Eq. ( |162| ) we see that the jet mass, mj(fi), depends on the renormalization scale. 
This dependence arises because the jet functions have an anomalous dimension. In the jet- 
mass scheme we induce additional /i-dependence in the cross-section that appears through 
mj(fi) in the variable s and through the Smj(fi) term in B±. This /i-dependence cancels 
out order by order in perturbation theory. To implement the jet mass in the factorization 
theorem we proceed as follows. We take the value of the jet mass at a certain reference 
scale, mj(fMo), as the parameter one would like to determine from fitting the cross-section to 
data. In terms of this parameter one determines the jet mass mj(fir) via Eq. (|163[) where 
fi = fir is the scale for which the jet function is to be determined perturbatively. In the jet 
functions B±(s, T t , fi r , 5m) we then use 5m = 5m j (fir) and the invariant mass variables s t ,t 
are determined via 

Mf T -m 2 j(fir) 

hi = ' , ■ (165) 

mj(fir) 

We emphasize again that it is crucial that the renormalization scale in 5m and in the explicit 
logs in the jet functions B± agree, in order to ensure the cancellation of the (9(Aq CD ) renor- 
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malon ambiguity. This is because asymptotically (for large order n) the terms causing the 
poorly behaved perturbative behavior in the pole mass scheme at (D(a™ +1 ) are proportional 
to [/ir a™ +1 (/ir)/3o ra!], and only cancel if the same renormalization scale is used. 

Relation to other mass schemes. It is useful to relate the jet-mass to the top MS mass, 
m t (fi). This facilitates using a top mass measurement from jets in other computations, such 
as electroweak precision tests. Typically one is interested in m t (m t ), since the renormaliza- 
tion group evolution in MS makes sense only above the mass of the particle. To relate the 
two mass schemes we take the measured mj(/i ) and use the solution to the jet-mass RGE 
equation in Eq. ( |163| ) to run it up to, lets say, [i = m t , obtaining mj(m t ). Now we use the 
relations to the scale independent pole mass at /i = m t : 

c(m t ) m t (m t ) = m p t oXe = mj(m t ) + Cj(m t , L m ) L m , (166) 

where c(/i) = 1 + CVa s (/z)/7r(l + | ln(/i/m t )) + . . . and cj is given to one-loop by Eq. (|162[) . 
Recall that the choice of L m determines a scheme, so rrij(n) also depends on this parameter. 
Expanding the relation in Eq. ( |166| ) to one-loop order we obtain a translation of the jet-mass 
to the MS scheme that is free of the (9(Aqcd) renormalon contained in the pole mass: 



/_x a s (m t )C F a s (m t )C F r fm t \ 3] 
m t (m t ) = mj(m t ) mj(m t ) H L m In — + - . (167) 

7T 7T " 



V Lm ' 2 



Note again that it is essential to strictly expand the series on the RHS of Eq. (167) and to use 
the strong coupling constant at the same scale everywhere to ensure the proper cancellation 
of the C(Aqcd) renormalon ambiguity |9T|. We also note that c(m t ) is known to three-loop 



order |)^, |93|, |94|, while cj(m t , T t ) is only known to one-loop at this time. Due to the small 
size of T t the one-loop contribution of cj(m t , T t ) causes only a shift of about ~ 250 MeV in 
the determination of the MS mass fnt(mt). Thus this correction may in many cases not be 
of critical concern when converting a top-mass determination from jets into an MS mass at 
one-loop order. However, we emphasize that mistaking a jet mass measurement as a pole 
mass value beyond the one-loop order can lead to a significant error in precision quantities 
that have a strong dependence on the top quark MS mass and which have been computed 
to high order in QCD. This is due to the (9(Aqcd) renormalon inherent to the pole mass 
definition. It is therefore an important task to determine the higher order contributions in 
the jet mass definition of Eq. (|162|). 



Another important class of top quark masses are the so-called threshold masses |80 
which can be determined to very high precision from a threshold-scan of the total top pair 
production cross section at a future e + e~ linear collider. Based on theoretical predictions at 



the next-to-next-to leading level in QCD [^, |8T], |95| and through dedicated experimental 



studies it is expected that a threshold top mass such as the lS-mass [59, B9L M can be 



determined with theoretical and experimental uncertainties at the level of 100 MeV 0, |97 



It is therefore useful to relate the jet mass to the top IS mass. To establish this relation 
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one should note that the IS mass is defined in the framework of nonrelativistic QCD and 
incorporates effects which are associated to soft ~ m t a s and ultrasoft ~ m t a 2 s scales. Since 
ultrasoft effects are not responsible for the nonrelativistic binding effects that define the 
IS mass definition and since the (9(Aqcd) renormalon contribution in the lS-pole mass 
relation are associated to the soft scale the relation between the jet and IS mass has to 
be determined for the soft scale /is ~ mtOt s . To obtain the relation one can use an approach 
similar to the one described above and first evolve the jet mass to /is using Eq. ( |163| ). It is 
then straightforward to relate the jet mass to the IS mass using the known results for the 1S- 
pole mass relation, see e.g. Refs. |81j, |98], for three-loop results (see also Refs. |100| , |101|| ) 
accounting also for summation of large logarithmic terms and Refs. [ |102| , 103 ] for four-loop 
fixed order expressions. At one-loop order the relation reads 



is 



8 



Oi s (fJ>s) C F mj(/iL 



. a s (/i s )C F 

H J^r, 



(168) 



Note that the same principles for treating the perturbative series discussed above for the MS- 
jet mass relation have to be applied here to ensure the proper cancellation of the C(Aqcd) 
renormalon contributions. In addition it is necessary to treat the terms in the perturbative 
series in the lS-pole mass relation in the so-called Upsilon expansion, where terms that are 
of order are formally treated of order a n s |89|, This is because the physical scale 
that governs this series is the inverse Bohr radius Cpm t a s (which is the analog of L m in 
Eq. ( |162| )). Note that the one-loop corrections from the jet-pole mass relation have a larger 
numerical impact in the one-loop relation of Eq. ( |168j ) than in Eq. ( |167| ) because the lS-pole 
mass corrections are an order of magnitude smaller than the corrections in the MS-pole 
mass relation. We will give a more detailed discussion on the higher order structure of 
Eqs. SJHD in Ref. 



VII. SOFT FUNCTION MODELS WITH PERTURBATIVE CORRECTIONS 



The soft function at a scale /i ~ /i\ is written as 



+oo 



Spart (£ + -t, £--£-, fl, 6) S mod (I + ,t) 



(169) 



This combines the partonic perturbative result for the soft function 5 ,part (given in Eqs. ( |113| ) 
for the hemisphere prescription), with a model hadronic function S mod satisfying the moment 
constraints in Eq. fl4ip. As explained in Ref. |HJ, this form encodes the features we require 
for an appropriate soft-function 5* for our analysis. In particular it works equally well for the 
peak region where the soft-function in non-perturbative, and for the tail region where the 
soft-function is perturbatively calculable at leading power. S in Eq. (|169|) has /i dependence 
consistent with its anomalous dimension and the MS scheme. And finally it should be totally 
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free from the (9(Aqcd) soft-function renormalon ambiguity identified in Ref. [30], which is 
also known to appear in event shapes for massless jets ||104|| . 

For the analyses in this work we will use the exponential model / ex p of Ref. [52], with the 
addition of a gap parameter A, so that 



'S'mod {£ j 



A) = / ( 



cxp \ 



fexp{£ j £ 



A,r-A), 

N JV(a,6) 



(170) 



a-i /-(£+) 2 -(£-) 2 -26£+r 
exp 



' A 2 V A 2 J ~r V A 2 

Here the normalization constant Af(a,b) is defined so that J d£ + d£~ S(£ + , £~) = 1. The 
parameter A ~ Aqcd sets the width of the hadronic function and hence the scale for £ ± and 
the soft radiation. The dimensionless parameter a controls how fast the function vanishes at 
the origin, and the dimensionless parameter b > — 1 controls the correlation of energy flow 
into the two hemispheres. Any b ^ implies cross-talk between the two hemispheres. 9 The 
gap parameter A enforces I > A and encodes the minimal hadronic energy deposit due to 
soft radiation. 

As explained in Ref. [[HJ, there is a renormalon in 5 , part (^ ± — £ ± ) that corresponds to 
an (9(Aqcd) ambiguity in the partonic threshold where I — £ = 0, and a corresponding 
ambiguity in the non-perturbative gap-parameter A. It can be removed by shifting to a 
renormalon free gap parameter A, using = A = A(/x) + 5(fi), 



s(£ + , r,n) 



oo J — oo 

+oo p+oo 

~+ 



dV 



d£~ S paTt (£ + -£ + ,£--r,fi)f cxp (£ + -A,£--A) (171) 
S part (£ + -£ + -5, r-t-5, fj) f exp (£ + -A, r - A) . 



OO 



Here 5 = $i * s a perturbative series with 5i ~ 0(a l s ) that defines the scheme for A. 
Expanding 5 , part (i' ± — £ ± — 5) in perturbation theory the <Vs remove the renormalon ambiguity 
from S'part order by order. Up to 0(a s ) this gives 



S pa ,rt(£ + ,£ T^i^i) — 'S'part (^ + 1^ ) + 



( s; art (£ + ,r, /U )-5 1 



d d 



^~J7T Impart (^ + > ^ ) 



,d£+ d£~J part 
where defining C l {£) = l//x \ 9{t) ln(£///) /(£///)] we have 

SpU^O = 5{£ + )5{£-) , S^{t,tr,ti = 5{£ + )Sl Tt {£- ,/i) + <KO^ art (£+, //) 



(172) 



'S'part (A A 4 ) 



C F a s (/j,) 
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A renormalon free scheme for the gap A can be defined |30| using a first moment of the soft 
function with upper cutoff La, similar to the jet-mass in Eq. (|157|) . This definition can be 



9 In Ref. |Q the values a = 2 and b — —0.4 were obtained from a fit to LEP data. The analysis used a 
different scheme for including perturbative corrections in the soft-function than the one advocated here. 
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FIG. 11: Soft function models based on Eq. ( |178[) with the hadronic model function S mod given 
in Eq. O for A = 0.55 GeV, A(/i = 1 GeV) = 0.1 GeV and (a, 6) = (2.5,-0.4) (left panel), 
(3.5,-0.2) (middle panel), and (a, b) = (2.5,-0.8) (right panel). The curves are tree-level (black 
solid line), 0(a s ) with <5i = (blue dashed lines) and 0{a s ) with a renormalon free gap (red light 
solid lines). The blue and red curves are shown for fi = 0.8,0.9, 1.0 GeV, with the higher curves 
correspond to lower values of the renormalization scale. 



written 



A di + [ A d£- t S paTt (£ + - 5, r - 6, /i) , 



and at 0(a s ) gives 



5i 



-2L, 



C F a s (fi) 



TV 
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(174) 
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Because A = A(/i) + 5(fi) is RG-invariant, this gives an anomalous dimension equation 

*A {p) = 2L.£^, (176) 
an n 

with a LL solution 

^ ra '^l (177) 



4CV 



Using Eq. (|172|) in (|170|) and integrating by parts we obtain a suitable soft-function for 
our NLL analysis 



S(t , r , n) = S mod (£ + , r , A) - 5 1 ( 4t + 4r ) <s i > 



A) 



(178) 
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oo J — oo 
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si Btt (£ + -£ + ,r-r,n)s mQi (?;r,A) 

£ part (£+-£ + ,/i,5i) S^itr-t,^) S mod (£ + J-,A) , 

where the modified one-dimensional partonic soft-function is 

S part (£, /i, tf x ) = 5(£) - 5x(/i) ) + ^ art (l, /i) . (179) 
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With Eq. (|170|) the result in Eq. (|178[) involves logarithmic terms In (£///) and ln(A//x), 
that arise from the convolution of the partonic plus functions with the smooth hadronic 
functions. In the peak region there is a possible tension between the convergence of the 
perturbative series in a s (n) for S part and the size of the ln(A//x) terms. In Ref. [BO] the log- 
series in the soft-function was analyzed and it was concluded that for scales \x ~ 1 GeV in 
the peak region this tension is not an issue. In the tail region the cross-section is dominated 
by ^ ± ~ sm/Q. Thus grows, and it becomes necessary to increase /i so that fi ~ sm/Q 
to avoid large logs from the In^/n) terms. 

To demonstrate the importance of the renormalon subtraction we show S(£ + ,£~,fi) in 
Fig. [O], plotted with £ = £ + = £~ for A = 0.55 GeV and three sets of the remaining param- 
eters (a, b, L A /A) = (2.5,-0.4,0.8), (3.5,-0.2,0.7), and (2.5,-0.8,0.8) respectively. The 
tree-level soft function is S = S mo( ± in Eq. ( |170|) where we take A = 100 MeV, and is shown 
by the black solid lines. The three dashed blue lines denote the 0(a s ) soft function obtained 
from Eqs. (|178| ) without the renormalon subtraction (Si = 0), and for fi = 0.8, 0.9, 1.0 GeV. 
The three light solid red lines denote the 0(a s ) soft function obtained from Eqs. ( |178j ) with 
a renormalon free gap using Si from Eq. (|175|) , and \i = 0.8, 0.9, 1.0 GeV. For both the blue 
and red curves we use A(/i) = 60, 82, 100 MeV for these three /i's respectively, which corre- 
sponds to implementing the LL running from Eq. ( |177| ). Compared to the tree level result, 



the blue dashed 0(a s ) curves show a significantly shifted location of their maximum and 
become negative for small values of £. The red light solid curves show that the renormalon 
subtraction stabilizes the peak location and removes the negative dip. These features are 
generic for any choice of a hadronic model parameters (a, b). 

VIII. NUMERICAL ANALYSIS UP TO NEXT-TO-LEADING LOG ORDER 

In Ref. we carried out a numerical analysis of the top-invariant mass distribution 
concentrating on nonperturbative effects caused by the soft function and on the dependence 
of the invariant mass distributions on the parameters used for the soft function model. 
The analysis was based on the tree-level results for the jet functions in Eq. (|28|), without 
summation of logarithms, and on the hemisphere soft function as given by the model of 



Ref. [22] which had been obtained from fits to event shapes in e + e~ annihilation. The 
soft-function caused a positive shift in the peak position of the invariant mass distribution, 
M tj f > mj, where the shift is parametrically ~ AqcdQ/^, and it was demonstrated that 
the peak shift and peak width grow linearly with Q/m. 

In this section we will extend the analysis to include radiative corrections and examine the 
perturbative convergence of predictions for the invariant mass distribution. This amounts 
to a full NLL analysis (i.e. one-loop matrix elements plus NLL summation of logarithms). 
Recall from Fig. [I] that there are four relevant scales for the log-summation, /j,q ~ Q, 
fi m ~ m, fir — s + T t + QA/m t and fi\ > A QCD + mT t /Q + sm/Q. For the peak region 



59 



we use /iA — 1 GeV. Our analysis is performed in several steps. After setting up the 
cross-section formula, we proceed in section |V111 A| to consider the summation of large logs 
for the perturbative corrections, and analyze the scale and scheme dependence. We show 
that the summation of logs between /ir — T t + QA/m t and fi\ ~ 1 GeV have a significant 
impact on stabilizing the cross-section. Then in section [VT l l~B| we convolute the perturbative 
corrections with the soft-function model and analyze the cross-section in the peak region. In 
section |V111 (J| we analyze the cross-section in the tail region, and plot combined peak and 
tail results. Finally in section |V111 D| we use our results to determine the thrust distribution 
at NLL order. 

For the numerical analysis it is convenient to write the invariant mass cross-section in 
the top jet-mass scheme in terms of dimension one invariant mass variables 



d 2 a 



^F(M t ,M h mj,-Q-), (180) 



dM t dM- t Tf V ' " " mj 

where the pref actor o"o is given in Eq. (^). Here mj is the jet-mass and the dimensionless 
function F is 

F (M t , Mt, mj, — ") = / V r P - ^, % - ^, M ) 5™ d (£ + , r , A( MA )) (181) 

30 ^ m J m J m J m J ' 

with 5' mod the hadronic model function given in Eq. ( |1 70|) . In the second line we shifted 
the integration variables to put all /i-dependent factors into P. In terms of M t and Mj the 
invariant mass variables s^j in Eq. ( |181| ) are 



M'f - mi Ml - mi 

s t = — , s { = — l - J - . 182 

mj mj 



All the perturbatively computable contributions in Eq. ( |181| ) are grouped into the dimen- 
sionless function 

P (s t , §t, // A ) = 4:M t MtT? H Q (Q, h q ) U Hq (Q, h q , /i m ) H m (mj, fx m ) U Hm (^-, fx m , /i A 

x G + (s t , — ,r t ,/i A N )G_fs i -, — ,r t ,/i A ) . (183) 
V mj / V mj / 

P also depends on Q, mj, and T t , but for simplicity we have not shown this dependence in 
its arguments. For the hard coefficients we used Eqs. (|55|) and (65) to write them in terms of 
the one-loop matching coefficients Hq^Q^q) and H m (m, /z m ) in Eqs. (|93|J132|) and the NLL 
evolution factors Uh q and Un m given by Eq. ([791). The functions G± in Eq. ( |183|) contain 
perturbative corrections that modify the shape of the cross-section. Using Eqs. (|46|j39| , p!78 
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and a few trivial changes of integration variables, these functions are 

( Q \ f +oc 

G±[s, — ,T t , Ua) = / ds' ds" d£' Ub(s — s' ua, ttr) 

x £>r=o _; S »_Q_ g ^ 6m \ § /iA; 5 ) F * . (184) 

V mj / 7r(s"^ + 1 1 ) 

This result depends on -B± =0 , the jet function for stable quarks in Eq. ( |142|) , and 5 ,part the 
modified partonic soft function of Eq. ( |1 79| ) . The form in Eq. ( |181| ) is derived from the 
factorization theorem given in Eq. ([46]), where the renormalization scales /ir and fi\ were 



distinguished. This leads to the presence of the evolution factor Ub in Eq. (|184|) , which is 
given at NLL in Eq. (170) ■ The functions G±, and hence all the ingredients in P, can be 
computed in perturbation theory, and analytic results for G± are given in Appendix [Ej. 

When quoting results at LL order we take Ub, Uh q , and Un m at LL order, and use tree- 
level results for £?£ =0 and S pavt , including 5m — 8\ — 0. The results quoted at NLL order 
use NLL-evolution for Ub, Uh q , and Ujj m . They also include the 0(a s ) results for matching 
coefficients and matrix elements, including B± =0 , S pax t, Sm, 6%, Hq(Q,hq) and H m {m, /i m ). 
These 0(a s ) terms have no-large logs, and in our numerical analysis we strictly drop all 
terms of 0{a 2 s ) or higher in the product of these matching and matrix element terms that 
appear in P. We also make use of the two- loop solution for the running coupling 



+ r ln (-) + tt ln I 1 + r a -M ln (-)l > (185) 



a s (/i) a s (/i ) 27r V/i / 47t/? 

with a s ([io = mz) = 0.118 as our reference value, and with (3q and (3\ from Eq. (p3|) . For 
the running above /x m we take n/ = 6, while for running below /i m we take n/ = 5 (hence 
neglecting the 6-quark threshold). 

Since there are many features of the cross-section formulae in Eqs. (p.80H184|) that we 
wish to explore, it is useful to have a default set of parameters to use at both LL and NLL 
order. When not otherwise specified, we use the following values for our analysis below. 
Our default Q/mj = 5, and the default renormalization scales are fig = 5 * 172 GeV, 
H m = 172 GeV, /ir = 5 GeV, and //a = 1 GeV. For results that involve the running jet mass 
we take as a reference value mj(/io = 2 GeV) = 172 GeV in the L m = 1 GeV scheme, and 
evolve to other scales using Eq. ( |163| ) for mj(fi). For results in the pole mass scheme we 
also use mf olc = 172 GeV. When the running gap is included we take as our reference value 
A(// = 1 GeV) = 100 MeV with L A = 0.44 GeV, and evolve to other scales using Eq. ([T77D 
for A(/i). We refer to this as the A-scheme. For results quoted without a renormalon free 
gap parameter we use A = 100 MeV, and refer to this as the A-scheme. Finally when 
studying perturbative aspects of the cross-section our default model for the soft function 



S mod (£ + , £~ , A) in Eq. is (A, a, b) = (0.55 GeV, 2.5, -0.4). Note that the dependence 
of S'mod on the model parameters A, a and b is not shown explicitly in its arguments. We 
will explore deviations from these default parameters on a case by case basis. 
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mfl+(s,r t ) S(£,£) P(M t ,M t ) 




FIG. 12: Renormalization group evolution of the jet function, mjB + (s,T t , fj,) (left panels), the 
diagonal soft function S(£,£,/j) (center panels), and the function P(M t ,M t ) in Eq. (|186| ) (right 
panels) with /i = = A*A- The top three panels show LL results, while the bottom three are NLL 
results. For the left and center panels the curves are tree- level (black solid), NLO at /i = 1 GeV (blue 
dot-dashed lines), and curves which evolve at LL or NLL order from [Aq = 1 to [i = 1.5, 4.0, 7.0 GeV 
(red, purple, magenta, with decreasing dash sizes respectively). The right panels show only these 
last three curves. 



A. Analysis of Perturbative Corrections and Log Summation 

We begin our analysis by studying the perturbative corrections contained in the function 

P [MuMl) , P (^-"S-^), . (186 ) 

V mj mj J 

This function P is convoluted with the soft-function model to give the cross-section, as shown 
in the second line of Eq. ( |181| ). All of the dependence on the scales /iq, /i m , fir, an d fiA 



cancels out in this function to the order we are working, and we can analyze the residual 
scale dependence of P to obtain an estimate of the remaining uncertainties from higher order 
corrections. 

In Fig. |TJ| we demonstrate the effect of summing logs for the functions mjB + (s,T t , fi) 
(left panels, see Eq. (|6^)), S(£,£,fi) (center panels, see Eq. (|59|)), and P(M t ,M t ) (right 
panels, see Eq. (|186|) ). For this plot we take fi = fir = fiA- The three top panels show LL 
results using the pole-mass and A-schemes. They demonstrate that increasing /i changes 
the shape of both B + and S. However, in the convolution that gives P these changes just 
reduce to a shift in the overall normalization, due to the consistency condition in Eq. ([78]). 
The three bottom panels show NLL results using the jet-mass and A-schemes. At LL and 
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FIG. 13: The an d fi\ scale dependence of the perturbative contributions, ¥(Mt,Mt). The 
top-panels show LL results, while bottom panels show NLL results. Central values are fir = 5 GeV 
and fi\ = 1 GeV. In the left panels the solid curves are for fir = 3.3, 5, 7.5 GeV (from top to bottom 
at the peak), while the blue-dashed line shows the result when fir = fi\ = lGeV. In the central 
panels the solid curves are for fi\ = 0.8, 1.0, 1.2 GeV (from bottom to top at the peak). The right 
panels are the same as the central panels, except that we also change fir so that fir / fi\ = Q / m = 5 
remains fixed. 

NLL order the peak of B + moves to the left as we increase fx (top-left panel and bottom-left 
panel). At both LL and NLL order the peak of S moves to the right for increasing fx (central 
panels). The right panels show that at LL and NLL we still have a strong residual scale 
dependence in P, and that there is still a peak shift at NLL. This occurs because we have 
taken fir = fi\ and not yet summed the large logs between fir and fi\, where fir/f^A — Q/m. 
These remaining large logs can be seen explicitly in the formula in Eq. (|E11|) . 

This situation is rectified in Fig. [13], where we show results for P(M t , M t ) with separated 
scales fir and fi\. Again the top three panels are LL, and bottom three are NLL. In the 
left two panels we vary fir about fir — 5 GeV holding fix — 1 GeV fixed. We use the 
range fir = 3.3-7.5 GeV to estimate the scale uncertainty because of the importance of 
not upsetting the fir/ Ha — Q/m relation too severely. (For contrast the blue-dashed curve 
in the lower- left panel shows the result for fir — fix — 1 GeV.) In the center panels we 
hold fir — 5 GeV and instead vary fi^ = 0.8-1.2 GeV (where variation below 0.8 GeV is 
not advisable since a s has grown to 0.45 at this scale). The dashed vertical line shows 
the input value of the short-distance mass, rrij(fi = 2 GeV) = 172 GeV, while the solid 
vertical line shows the shift due to the gap, Mj = [m 2 j(fio) + QA] 1 ^ 2 ~ 172.25 GeV, with 
A(l GeV) = 0.1 GeV. Finally in the two right panels we show how the fi- variation of the 
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FIG. 14: P(Mf, Mf) versus Mf. Left panel shows the fj, m dependence at LL (bottom curves) and 
NLL (top curves) taking \i m = 86, 172, 344 GeV. Central panels shows the \iq dependence at LL 
(bottom curves) and NLL (top curves) taking fiQ = 430, 860, 1720 GeV. The right panel shows the 
effect of using renormalon free gap and mass parameters, where the red solid curve includes both. 
The purple dashed curve turns off the renormalon subtractions for the mass (thus using the pole 
mass scheme), and the blue dot-dashed curve turns off the renormalon subtraction for the gap. 



center panels is reduced if we vary //a in the same range, but simultaneously change fir 
so that /ir//iA — Q/vn is fixed. Because of this sizeable correlation the overall uncertainty 
is smaller than from naively summing the uncertainties from the left and center panels in 
quadrature. We believe that an uncertainty in the shape that is of order the /ir variation 
shown in the left panels gives a reasonable error estimate. From comparing the percent 
change at LL and the percent change at NLL we see that there is a reduction to the \i- 
variation in all cases, particularly in the cross-section above the peak. 

In Fig. |T3| we show results for the fi m and [Iq scale dependence of P(M t ,M t ) (left and 
central panels respectively). We increase and decrease fi m and [Lq by a factor of two, and 
both variations exhibit rather small scale uncertainty. Here we show LL and NLL as the 
bottom and top three curves inside each panel. As is often the case in jet-physics, we note 
that there is a sizeable change to the normalization of the cross-section in going from LL 
to NLL order. The vertical lines are the same as in Fig. O. In the right most panel of 



Fig. [b| we show the effect of using the renormalon free jet mass and renormalon free gap 
parameter A (solid red curve) in contrast to turning off the renormalon subtraction for the 
mass, i.e. when using the top quark pole scheme (purple dashed curve), and turning off the 
renormalon subtraction for the gap (blue dot-dashed curve). Even at 0(a s ) the importance 
of having a renormalon free soft function with A is clearly visible. 



B. Cross Section in the Peak Region 



Having examined the scale dependence of the perturbative corrections we now 
turn to the convolution with the soft-function that gives the normalized cross-section 
F(M t , Mt,mj,Q/mj) in Eq. ( |181|) . For most of our plots we keep the soft function model 
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FIG. 15: F(Mt,Mf), the differential cross-section in units of 0o/rf , versus M t and Mj. The result 
is shown at NLL order. 




FIG. 16: Normalized peak cross-section, F(M t , M t ) versus M t . The dashed curves have /xr = 5 GeV, 
and the solid curves have /ir = 3.3, 7.5 GeV. The left panel shows results at LL (lower purple 
curves) and NLL (upper red curves) with the jet and A schemes. The center panel shows results 
in the jet-mass scheme (red) versus the pole-mass scheme (blue), where in both cases we use the 
A scheme. The right panel shows results in the A(/z) scheme for the gap parameter (red) versus 
the A scheme (magenta), where in both cases we use the jet-mass scheme. 



fixed, having in mind that it can be extracted from LEP data. In Fig. [15| we show F at NLL 
for our default parameter set as a function of the two invariant mass variables M t and Mj. 
The underlying short-distance quark mass is mj(/i = 2 GeV) = 172 GeV, and the peak of 
the cross-section occurs for M t and Mi values which are ~ 2.4 GeV larger. This peak shift 
occurs due to the presence of the low energy radiation described by the soft function as dis- 
cussed in Ref. 0. At LO the shift is in the positive direction to Mf esk ~ m J + QS^\/ (2m j), 
where here S 1 ' 1,0 ' = J d£ + d£~ £ + S mo( i(£ + , £~) ~ Aq C d is the first moment of the underlying 
soft-function model . As described below, this linear behavior with Q/m persists at NLL 
order, although the slope is no longer simply Sj^j- Above the peak one sees in Fig. [L5| the 
perturbative tails from gluon radiation, and that the tails are largest if we fix one of M t or 
Mi at the peak. 
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FIG. 17: a) Peak position of the single differential distribution da / dMtdMi(Mt, Mt) as a function 
of ^r- Red curves use the jet-mass and A scheme and blue curves use the pole-mass and A 
scheme. Dashed curves are LL order, and solid curves are NLL order. The bars to the left 
show the size of the scale variation and left-to-right correspond to the curves from bottom-to- 
top, b) The solid curves show the peak position versus Q/m for six different models, which from 
top to bottom are (a, b) = (3.5,-0.8) (purple), (a, b) = (2.5,-0.8) (blue), (a, b) = (3.5,-0.4), 
(magenta), (a, b) = (2.5,-0.4) (red), (a, b) = (3.5,0.4) (yellow), (a, b) = (2.5,0.4) (green). The 
solid curves show a linear fit using the values at Q/m = 4 and 5. Extrapolated to Q/m = any line 
converges on the underlying short-distance mass, independent of the soft-radiation model, yielding 
mtOu = 5 GeV) = 171.9 ± 0.1 GeV. 



In order to analyze the parameter dependence of the cross-section we will now consider 
the diagonal F(M t , M t , mj, Q/mj), which we simply referred to as F(M t , M t ) in the analysis 
that follows. In Fig. [Tj| in the left panel, we show LL curves (bottom three lines) and NLL 
curves (top three lines) using fir = 3.3, 5.0, 7.5 GeV in the jet-mass and A-scheme. We find 
that the peak of the cross-section is very stable to the variation of fir, and changes very little 
from LL to NLL order. As explained above, by far the dominant contribution of the shift of 
the peak away from the input short-distance jet-mass is due to the underlying soft-function, 
shown here by the difference between the dashed and solid lines. In the central panel we 
show again the NLL order cross sections in the jet-mass and A-scheme (red curves) and 
compare it to the NLL predictions in the pole-mass scheme for the same three fir values 
(blue curves). The results show that in the pole- mass scheme there is more variation of the 
peak position than in the jet-mass scheme. Finally in the right panel we show variations of 
the cross-section in comparing the renormalon free A-scheme (red curves) and the gap with 
a renormalon ambiguity in the A-scheme (magenta curves). This figure demonstrates that 
the effect of the switching to a renormalon free gap-scheme is larger than the residual fir 
dependence at NLL order. 

The fir dependence of the peak position is shown more explicitly in the left panel of 
Fig. [17]. In the pole-scheme (blue curves) we see that there is very little change to the fir 
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FIG. 18: Perturbative contributions, P(Mt, Mt) in the tail region at LL (dashed curves) and NLL 
(solid curves). In the left panel we take fi^ = arid plot three curves with fir = {0.5, 1.0, 1.5}iip. 
In the center panel we take fir = Mr and show three curves with fi\ = {0.5, 1.0, 1.5}/i A . In the 
right panel we show fir = {0.5, 1.0, 1.5}/u r with fi\ = fir m j(2 GeV)/Q- Here fi r and fi\ are given 



in Eq. (|i87|) . 



dependence in going from LL (dashed blue curve) to NLL (solid blue curve). In contrast in 
the jet-mass scheme (red curves) the fir dependence is already smaller at LL order (dashed 
red curve), and is significantly reduced by the NLL results (solid red curve). In the right 



panel of Fig. [17] we plot as black ticks results for the peak position versus Q/m using six 
different models for the soft function with fir = 5 GeV. It is clearly visible that the peak 
is shifted in a linear fashion with Q/mj, with a slope that is model dependent. For each 
model the solid lines show a fit to the the peak position for fir = 3.3 GeV and fir = 7.5 GeV 
with a solid band to show the uncertainty from this /i-dependence. The fits are done using 
the points at Q/m = 4 and 5. Extrapolating back to Q/m = removes the dependence on 
the soft-radiation, and we see that for any soft function model the intercept determines the 
short- distance mass parameter mj{fi = 5 GeV) = 171.9 GeV. From the spread of the curves 
we have ~ 0.13 GeV theoretical uncertainty in this determination of the short-distance mass. 
This provides a method for determining the short- distance mass even if the soft-function is 
unknown. In order to maintain the perturbative stability of the relation of this intercept 
with the top-mass it is important to use the jet-mass scheme. 



C. Cross Section in the Tail Region 



The tail region of the cross-section is characterized by invariant masses where s ^> T. In 
this region we are varying M t and hence s over a large range, and it becomes necessary to 
scale fir ~ s and fi\ ~ smj /Q to avoid having large logarithms that spoil the perturbative 
expansion. We therefore define reference scales for the tail region, 
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FIG. 19: Cross-section plotted over both the peak and tail regions at LL order (dashed) and NLL 
order (solid). In the left panel we take ha = Ha and plot three curves with ht = {0.5, 1.0, 1.5}/i r . 
In the center panel we take ht = Hr and show three curves with ha = {0.5, 1.0, 1.5}/t A . In the 
right panel we show ht = {0.5, 1.0, 1.5}/i r with ha = ht mj(2GeV)/Q. 



where mj(2GeV) = 172 GeV, and we study the scale dependence by varying /ir and fi\ 
about these results. In Fig. pl| we show the perturbative function P from Eq. ( |186| ) at LL 
order (dashed curves) and NLL order (solid curves). The left panel varies fir by 50% about 
/i° holding fi\ = fi A fixed, while the central panel varies /ia holding /ir = /^r fixed. In 
contrast to the peak region we now plot the cross-section over a log-scale. Note that the LL 
results exhibit larger uncertainty in this tail region, which is again substantially improved 
by the NLL results. In the right most panel we vary fi\ as in the central panel, but now take 
/W^a = Q/m as fixed. Just as in the peak region this choice substantially reduces the scale 
uncertainty, indicating once again that simply adding the individual variations of /ir and 
fi\ very likely overestimate the size of higher order perturbative corrections. Finally since s 
increases in the tail region, the uncertainty from the power expansion also increases as we 
require v ■ k/m = s/(2m) ~ (M t — mj)/mj <C 1. At M t = 185 GeV this is an expansion in 
1/12 and by M t = 200 GeV its an expansion in 1/6, both of which are larger than the ratio 
~ 1/100 that we have in the peak region. 

In Fig. [19] we convolute P with the soft-function model as in Eq. (|181[) , and plot the 
normalized cross-section F over both the peak and tail regions. The three panels show 
the same /z-variations as Fig. 18. These plots show one of the attractive features of our 
treatment of the soft-function. In the tail region there are perturbative corrections in the 
soft-function that are important for determining the cross-section, and there is also an 
important power correction due to the first moment of the model function S mo ^- Both of 
these are included in our analysis by using the full S(£ + ,£~,fi) function from Eq. ( |169| ) 
as explained in Ref. [M. In the peak region these terms naturally interpolate into a full 
model soft-function in a consistent manner. Once again we see from the third panel of 
Fig. [TP] that taking /ir/AU = Q/mj leads to quite small /i-dependence for the entire cross- 
section. Finally in Fig. ^Tj we show the effect that variation of the soft-function model has 
on the cross-section in the tail region. The effect of the soft function becomes larger as 
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FIG. 20: Variation of the tail cross-section using different soft function models at LL (left 
panel) and NLL (right panel). The curves are (a, b) = (3.5,-0.8) (purple), (a, b) = (2.5,-0.8) 
(blue), (a, b) = (3.5,-0.4), (magenta), (a, b) = (2.5,-0.4) (red), (a, b) = (3.5,0.4) (black), 
(a, b) = (2.5,0.4) (green). 



we get closer to the peak region, as expected. Since the cross-section has already dropped 
by two-orders of magnitude by M t = 200 GeV we have not bothered to analyze it in the 
ultra-tail region, M t — mj ~ mj, where it is further suppressed by several more decades. 
However, in appendix |F| we do give formulas for the cross-section in the ultra-tail region, 
which are analogous to the ones used for our analysis of the peak and tail cross-sections. 
These formulae could be useful as a means of estimating top-quark backgrounds from ti 
events for other processes in the ultra-tail region. 



D. Thrust 



Starting from the two-dimensional distribution, d a/dM 2 dM 2 in Eq. (|180|) it is straight- 
forward to derive results for other event shape variables for massive particles. For example, 
for the thrust T, we have 1 - T = r = (M t 2 + Mf ) /Q 2 , so 

1 da f°°, nr2 r° 2 j M 2 + Mf\ 1 (fa /100 , 
dM? I dMl 5[t - * ) -- 7T7¥7T7 t (188) 



a dT Jo Jo 1 V Q 2 ) a dM t 2 dM i 



V m J J 

The perturbative contributions are grouped into the dimensionless function Pt which is a 
projection of our function P, 

d (<* \ m jQ 2 v fs + s d s-s d \ 

Pt(*, = Jjs d P fi\j . (189) 

Here under the s d integral M 2 = m 2 + mj(s + s,i)/2 and M| = m 2 + mj(s — s^)/2. An 
analytic formula for Pt is derived in Appendix 0. The appropriate soft-function for thrust, 
S'Sod™^) in Eq. ( |188| ), is also simply a projection of the model for the hemisphere soft 
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FIG. 21: Thrust distribution, da/dT in units of ctq, plotted versus 1 — T at LL (dashed 
curves) and NLL (solid curves). In the left panel we take fi\ = fJl and plot three curves 
with fir = {0.5, 1.0, 1.5}^/p. In the center panel we take fir = Mr anc ^ show three curves 
with /iA = {0.7, 1.0, 1.5}/i A . In the right panel we show ^a = {0.7, 1.0, 1.5}/i A with fir = 
fi A [0.8 Q/mj (2 GeV)]. 

function, S mo d(^ + , &~, A), where 

/>oo 

^moT(^o)= / di + dr5(i-£ + -r)s mod (e + ,e-,o)- (iso) 

For the exponential model in Eq. ( |1 70| ) this projection gives 



where {a, b, A} are the model parameters and J\f(a, 6) is the same normalization constant as 
in Eq. fll70D . 

In Fig. |2l] we plot the thrust distribution at LL order (dashed curves) and NLL order 
(solid curves) for events which were initiated by the massive unstable top-quarks in e + e~ 
collisions. Since the plot includes values in the tail region we use the reference scales 



/ Q 4 / 2m 2 A2 2 o 1 fi T mj 

^ = V4^V " ~W) +( V) ' ^QA-Q-' (192) 

where r = 1 — T. Taking fir — fi r and fi\ ~ /x A ensures that the logs involving these 
parameters do not grow substantially over the region plotted. Our choice for fj,\ here is 
slightly larger than the ones used earlier. This is because of the effective doubling of the 
anomalous dimensions for the thrust cross-section (see Appendix |G|), which necessitates 
using slightly larger values for /z A to avoid the region where large values for a s cause a break 
down in perturbation theory. 

The threshold for thrust for two- massive particles is given by 1 — T = 2m 2 / Q 2 and is 
shown by the vertical dashed lines in Fig. [21]. Just as for the invariant mass-distribution, 
there is a peak in the thrust cross-section and it is shifted above the massive particle threshold 
due to soft-radiation effects by an amount ~ 2Aq CD /Q. The analog of this for massless jets 
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is a peak in the thrust distribution at values 1 — T ~ 2Aqcd/Q (see for example |22H), which 
is a shift above the massless dijet threshold at 1 — T = 0. The three panels in Fig. [H] show 
the //-dependence of our NLL results, varying /ir in the left panel, //a in the center panel, 
and H\ with /i r //iA = 4 fixed in the right panel (since here Q/mj = 5). Again we see that 
there is very small //-dependence when /ir and fi\ are varied in a correlated fashion. We 
believe the left panel gives a reasonable estimate of the perturbative uncertainties in the 
shape of the thrust distribution. An analysis of the thrust-distribution peak for different 
values of Q could also be used to extract the short-distance top-mass parameter. 

IX. CONCLUSION 

Precise measurements of the top quark mass m t belong to the most important standard 
measurements carried out at the Tevatron and the LHC The most sensitive method relies 
on the reconstruction of the top quark invariant mass distribution through measurements of 
the energies and momenta of jets from the top decay. While considerable work has and is 
being invested to control experimental systematic effects, very little theoretical work exists 
which studies both perturbative and nonperturbative QCD aspects of the resulting invariant 
mass distribution. Also, to our knowledge, there has been no theoretical work on how the 
shape and the resonance mass of this distribution are related to a short- distance top mass 
parameter in the QCD Lagrangian. 

In Ref. @ we derived the factorization theorem in Eq. (0) which describes the simpler 
environment of e + e~ collisions. It predicts the double differential invariant mass distribution 
cPa/dMtdM-j: in the resonance region for the large cm. energies Q 3> m t , where M t $ are the 
total invariant masses of all particles in the two hemispheres determined with respect to the 
event thrust axis. The factorization represents the leading order result in a power expansion 
in m/Q and T/m, and these corrections are indicated in Eq. (|^). Here T is the width of 
the invariant mass distribution, which is larger than the underlying total width of the top- 
quarks T t . The derivation was based on the hierarchy Q ^> m t ^> r,Aq CD , where Aqcd 
is the hadronization scale and uses the effective theories Soft-Collinear-Effective-Theory 
(SCET) and Heavy- Quark- Effective- Theory to achieve a separation of different physical 
effects associated to Q, m t , T t and M t f and Aqcd- For the systematic inclusion of mass 
and width effects the use of both effective theories was crucial. The factorization theorem 
separates perturbative from non-perturbative effects and represents the leading order term 
in the power expansion, but is valid to all orders in the expansion in a s . 

In this paper we extended the presentation given in Ref. and presented detailed com- 
putations of the different pieces entering the factorization theorem in the peak region at NLL 
order. We also presented NLL predictions for the tail of the invariant mass distribution, 
where M t) j are above the resonance peak. The double invariant hemisphere mass distribu- 
tion is itself an event-shape distribution that can be related to other event-shape variables 
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such as thrust or jet masses in a straightforward way. The factorization formula consists of 
several functions that can be computed perturbatively order-by-order in a s , including hard 
coefficients for the scales Q and m, and two jet functions for the top and antitop quarks 
which depend strongly on the top quark Lagrangian mass. It also involves a nonperturbative 
soft function that describes the momentum distribution of soft final state radiation. Using 
alternative invariant mass prescriptions, for which the soft particles are assigned differently 
to M t and Mj, the same factorization formula applies, but with a different soft function. In 
the tail region the soft function also contains perturbatively calculable corrections. 

Our analysis uses effective theory techniques. In particular we calculated Wilson coeffi- 
cients that arise from matching QCD onto SCET, and from matching SCET onto bHQET 
to NLO. In addition we calculated the NLL running: first between the scales Q and m using 
anomalous dimensions of operators in SCET, and then between the scales m and V and 
between V and A using anomalous dimensions of operators in bHQET. The perturbative 
corrections, including resummation, are given by simple analytic functions, and our strat- 
egy for computing these functions naturally generalizes for future analytic computations at 
NNLL order. One important result of our analysis is that in the peak region the running 
between the scales Q, m, and T is local, it only changes the overall normalization and not 
the shape of the invariant mass distribution. Thus only large logs between the scales T and 
A can shift the shape of the distribution. This is encoded in consistency equations for the 
renormalizat ion-group evolution functions entering the factorization theorem. 

The observable invariant mass distribution is obtained from a convolution of the jet 
functions with the soft function. Through this convolution the energy of the peak and 
the width of the observed distribution are dependent on the center of mass energy Q. In 
particular, non-perturbative effects described by the soft function shift the resonance peak 
position towards larger masses, and broaden the distribution. In Ref. [@ we demonstrated 
that the soft function for the hemisphere mass prescription can be determined from event 
shape distributions for massless dijet events, and that the dependence on the mass scheme 
is controlled through the perturbative expansion of the jet function. This allows in principle 
for a top quark mass determination that is free of hadronization uncertainties. Even if the 
soft-function is not known from measurements of massless jets, one can still extract the 
short- distance mass parameter using an analysis like the one described in Fig. |17|. We have 
demonstrated that these statements remain true in the presence of perturbative corrections 
and with the summation of large logs. We also introduced the so-called jet-mass scheme to 
define the Lagrangian top mass parameter. This is a top quark short-distance mass scheme 
that is particularly suited to mass determinations related to the resonance peak position, 
since it makes the peak position stable to the inclusion of perturbative corrections. One-loop 
relations of this jet-mass to the MS-mass and lS-mass schemes were given in section |VT|, as 
well as the LL evolution formula for the jet-mass. For the construction of the soft function 
at NLL order we used results from Ref. [3D] where a soft model function is convoluted 
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with the soft function contributions determined from fixed order perturbation theory. This 
soft function works equally well in the peak and tail regions. To avoid large higher order 
corrections it was necessary to introduce a gap parameter in the soft-function that accounts 
for the fact that there is a minimal hadronic energy for the soft-radiation between the jets. 
This parameter allows us to make the perturbative corrections in the soft function free from 
an (9(Aqcd) renormalon ambiguity. 

In our numerical analysis we analyzed the hemisphere mass distribution on the peak 
and away from the peak, and the thrust distribution. We demonstrated that NLL order 
corrections are important and need to be accounted for to make viable predictions. We also 
showed that it is important to sum the large logs between fir — T and /xa ~ A to avoid 
sizeable scale uncertainties. We also studied the impact of the jet mass scheme, and showed 
that it improves predictions for the resonance peak position in comparison to the pole mass 
scheme. For mass measurements this result implies that the jet mass can be determined 
from mass reconstruction more accurately than the pole mass which is known to suffer 
from renormalon ambiguities. In our NLL analysis we demonstrated that the perturbative 
corrections associated to the gap parameter improve significantly the perturbative behavior 
of predictions in the peak and in the tail region. This result implies that soft functions that 



account for a gap as proposed in Ref. [3D] are crucial for precise measurements of the top 
quark mass and the model parameters from experimental data. Finally we also presented 
a NLL prediction for the thrust distribution for top pair production in the region of large 
thrust. To our knowledge this NLL result presents for the first time a full resummed event- 
shape distribution for massive quarks. The thrust distribution has a strong dependence 
on the mass and can serve as an alternative way of measuring heavy quark masses. Our 
numerical analysis can be extended to make predictions for bottom quark production by 
taking the limit T t — ► 0. It would also be interesting to study the a s m t /Q and m 2 /Q 2 power 
corrections, which can be accomplished in our effective field theory setup. 

Through our detailed calculations of the jet invariant mass distributions and their relation 
to the top quark mass, we have demonstrated the viability of extracting the top mass 
with high precision at a future linear collider such as the ILC In principle, a precision 
of better than Aqcd can be achieved since there is a clear relation between the top mass 
Lagrangian parameter and the physically observed jet invariant mass distribution. In the 
future we intend to extend the work presented here and in Ref. [[| to the study of top-mass 
reconstruction at the LHC 
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APPENDIX A: SUMMARY OF FEYNMAN DIAGRAMS IN SCET 



In this appendix we list results for the individual Feynman diagrams used in the body of 
the paper in section [TV[ 

QCD graphs. For the QCD current at one-loop we have vertex and wavefunction graphs. 
We use free quark external states with an offshellness IR regulator, A 2 = p 2 — m 2 = p 2 — m 2 
where A 2 -C m 2 . We use dimensional regularization with d = 4 — 2e for the UV divergences. 
The graphs in Fig. are 
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SCET vertex graphs. For the SCET current we have collinear vertex and wavefunction 
graphs and a soft vertex graph. We use dimensional regularization for the UV divergences 



and the offshellness IR regulators A 2 = p 2 — m 2 and A 2 



p 2 
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To compare IR 



divergences to the full theory results we should set A 2 = A 2 and expand in A 2 <C m 2 . For 
later convenience we will first quote results prior to making this expansion. The massive 
collinear quark Lagrangian is given by ||105| , |106 
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Using the corresponding massive SCET Feynman rules the graphs in Fig. |5| are given by 
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where 
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(A4) 



and where all terms in the denominator with square brackets are defined with the +i0 
prescription. To evaluate Fig. [5]a,b we included the zero-bin minimal subtraction [|55| to 
avoid double counting the region encoded in Fig. [5]c. For and we have a singularity 
for n ■ (k + p) — > with fixed kj_, but it cancels in the sum of the two diagrams: 
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Setting n ■ p = n ■ p = Q and computing the integrals we find 
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Note that the soft graph, V| c , is independent of m. Expanding V| a and V|d + V|J= for A 2 <C m 2 
gives 
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From Eq. ( |A2j ) with Eq. (pT[ ) we have a mass counterterm, —i^,/(2n-p)2m5 m . In the pole 
mass scheme 
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which exactly cancels the entire m 2 term in the self energy graphs 
mass scheme we have S m = (5 m ) pole + 5m, and 
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Using the MS subtraction, the wavefunction renormalization removes the remaining 1/e 
divergence in Eq. ( |A8|) , and leaves a finite correction to the residue of the collinear quark 
propagator, iR^(^/2)n-p/(p 2 —m 2 ), with 
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This is identical to results from the QCD self energy, Vfck) in Eq. (|Al|), namely = R^ 



and Z 



SCET jet function. The computation of the stable massive SCET jet function is given 
by the imaginary part of the graphs in Fig. ||]. The integrals are identical to those for the 
vertex graphs above, except that now there is true external momentum, r^, that is routed 
through the diagram, and the invariant mass, s, takes the place of the offshellness, therefore 



these diagrams do not require an IR regulator. Here s 



m 



Qr£ — m 2 is defined for 



convenience. The diagrams in Fig. |] contribute to the SCET jet functions. Taking the spin 
and color trace, but not yet the imaginary part, the graphs give 
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where 1/e terms are UV divergences and all s factors are s + iO. An internal Z^ counterterm 
is not needed since these factors cancel between the propagator and vertices. The sum of 
the terms in Eq. ( |A10| ) is quoted as Eq. fl99|) in the text. Expanding Eq. (|A10|) to leading 
order in s/m 2 <C 1 gives 
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and the sum of these five terms gives Eq. (fUXj). 

Soft function graphs. Next we summarize the computation of the hemisphere soft function 
at one-loop given by the graphs in Fig. ^ We use dimensional regularization for both UV 
and IR divergences. For Fig. [7|a,b we have a loop integral, and for Fig. 0c, d a phase space 
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integral: 
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Here and S|j]b are scaleless and convert IR divergences in SfLd into UV divergences (see 
for instance Ref. [J7!J where this is worked out explicitly in several cases). To integrate 3^ t d 
we convert d d ~ 1 q = 7r 1_e /[2r(l— e)] dq + dq~(q + q~)~ e (q + +q~)9(q + )9(q~). Evaluating the sum 
of diagrams we find 
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This result is used in Eq. (|108|) of the text. 
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APPENDIX B: SUMMARY OF FEYNMAN DIAGRAMS IN BHQET 

In this appendix we list results for the individual Feynman diagrams used in the body 
of the paper in section [V|. The velocity four vector v± and the momentum fluctuation four 
vector k± for the top and antitop respectively are given by 

«?=f£.2.oA * f ~r("«iY (Bi) 



Q m J ~ \Q m 



m Q J ~ \m Q 

The Feynman rules for boosted HQET are summarized in Fig. [22| 

bHQET vertex graphs. First we give the results for the bHQET vertex graphs of Fig. || 
These graphs are defined in terms of the integrals I^ a ,6,c,d as Fig.||K = h v+ T^h v _) for 
k = a,b,d and Fig.|[r = for r = c, e, /. The one loop integrals are given by 
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soft 



= igT A n^ 



ucollinear 




mv+ + k 
> 



v + -£ + ,T r /2 



FIG. 22: bHQET Feynman rules for the top quarks annihilated by h v+ : soft gluon coupling to the 
top quark, ucollinear gluon coupling to top quark, Wilson line ucollinear gluon coupling to the top 
quark, top propagator. Results for the antitop quarks annihilated by h v _ are obtained by taking 
t> + — > V-, T A — ► = — (T^) T , and n «-> n. 



where all the external legs are kept off shell by (v± ■ r±) ^ to regulate infrared divergences. 
As before, throughout this section all factors in the denominator with square brackets are 
denned with a +i0 prescription. Just as in SCET, the integrals are defined with a zero- 
bin subtraction [53] in order to avoid double counting the region encoded in the soft loop 
in Fig. [|c. The loop momentum appearing in V^^^ has ultracollinear scaling as displayd 
in Eq. (]B1|). On the other hand, the loop momentum k in has a homogenous soft 
scaling corresponding to the exchange of a soft cross-talk gluon between the top and antitop 
sectors. Given homogenous scaling of the soft loop momentum and the collinear scaling of 
the velocity labels in Eq. (|B1|), at leading order we have the relations v + ■ k = \{n- v + ){n ■ k) 
and V--k = ]?(n-vJ)(n- k) which have been used in VEL Using v± -v± 
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From VjgJ; for the wave function graph we get the MS wave function renormalization Zh and 
the corresponding residue Rh as 
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These results are used to obtain Eq. (|128 ) in the text. 
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bHQET jet function. Next consider the graphs for the stable bHQET jet functions in 
Fig. |9|. Prior to taking the imaginary part they are 
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From Eqs. (|B 



and (p3| ) we get 
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where s = s + iO. Adding together 
body of the paper. 

Zero-bin subtraction for the bHQET jet function. The 1/e singularities in Eq. ( p6|) are 
UV divergences. This is ensured by the zero-bin subtraction [^] for the graphs B^"® that 
are needed to avoid double- counting with infrared regions already accounted for by the soft 
function (or the contributions in the soft loop in V|c). Using dimensional regularization 
to regularize UV and IR divergences in Eqs. ( p2|) this may not be obvious since in this 
case the zero-bin subtraction is associated to scaleless integrals. To illustrate the role of the 
zero-bin subtraction more explicitly we reconsider the calculation of the anti-top jet function 
B^" with a different IR regulator. From the definition of B+ =0 in the effective theory this 
computation involves two terms, the naive loop integrand and a term induced by the 

zero-bin subtraction on propagators -^^(o)- ^ e use an explicit p-term to regulate the soft 
and collinear IR divergences for n ■ k = k + — » 0. For the naive part of the result we have 
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Tim s An 

For the graph -B^ the zero-bin subtraction is obtained from the fact that the collinear 
propagators act as distributions and induce a subtraction from the limit where k^ ~ TA 2 . 
Since this subtraction is obtained from the Feynman rules of the same diagram as Eq. (B7), 
it must have the same p-regulator for the k + — > limit. For this subtraction part we obtain 
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The p term in Eqs. ( |B7|) and ( [B8|) regulates soft and collinear divergences, and these two 
results can be computed explicitly in terms of this regulator. However it drops out in the 
difference of Eqs. (|R7|) and (|B"8|) which is IR finite: 
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For k + — > the term in the brackets is linear in k + and the p-term can be dropped since the 
expression does not contain any IR-divergences. Evaluating Eq. ( |B9p we recover the result 
shown in Eq. Thus all 1/e singularities are indeed UV divergences. 



APPENDIX C: PLUS FUNCTION AND IMAGINARY PART IDENTITIES 



The plus function with an arbitrary exponent 1 + uj with uj < 1 is defined by 
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A more general definition for the distribution [6(x)/x 1+UJ ] + can be defined, by first integrat- 
ing 8(x)/x 1+UJ with a test function for values of uj where the integrals converge, and then 
analytically continuing the result to other values of uj. Expanding Eq. ( |Cl| ) for small to gives 
the definition for the log n (x)/x plus-functions for n > 0: 
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The following rescaling identity for a dimensionless constant k is also quite useful: 
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For example, for integrations over a finite range Eq. ( |C3|) allows us to rescale the +- function 
to act only within the interval [0, 1], where standard identites such as J Q dx'g(x') [9(x')/x'] + = 

dx'[g{x') —g(0)]/x' for a given test function g can then be applied. This is somewhat sim- 
pler than the corresponding general relation L dx'g(x')[6(x')/x'} + = L dx'[g(x')—g(0)]/x'+ 
g(0) ln(/3). Relation (|C^ ) has also been used to verify the multiplicative form of the consis- 
tency conditions in Eqs. ( |121| ) and (|150|) . 

In computing the SCET and bHQET jet functions one has to take the imaginary part of 
the forward scattering graphs. For this the following result for a dimensionless variable x is 
quite useful 
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where [p] on the sum is the greatest integer not exceeding p, sometimes also called the Gauss 
bracket of p. For the first few orders this gives 
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To compute the massive SCET jet function in Eq. ( |101| ) the following identities were also 
used [with s = s + iO = xnf + zO]: 
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APPENDIX D: GENERAL RGE WITH PLUS AND DELTA FUNCTIONS 



In this appendix we solve the general anomalous dimension equation 
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where t and t' are variables of mass-dimension j, and 7f(£ — t') involves a +- function and 
5-function. To motivate the general form for 7 we consider a generic 1-loop amplitude which 
has the form 
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At one-loop the imaginary part of y4 bare is renormalized by the ultraviolet Z-factor 
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where the numerical coefficients To and 70 are the first terms in the perturbative series for 
the anomalous dimensions 
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At any order in a s the anomalous dimension to be used in Eq. ( pi|) is 
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and depends on the dimension j, t — t', and /x. For convenience we introduced in Eq. ( p5| ) 
the mass scale k > so that the plus function has the dimensionless variables tj^ and 
t' j k? . But note that jf(t — t', /i) is independent of the choice of k. Here T[a s ] and j[a s ] are 
perturbative series in a s (fi) which start with a linear term as shown in Eq. (P4|) . 

To solve Eq. (P5|) we use the Fourier transform method of Ref. p3j 107 ], including the 



improvements of Ref. ||108|| which gives formulas that apply to all orders in perturbation 
theory. Our computation is a simple generalization of these solutions to mass- dimension 
j variables, which as we will see is key to understanding how the renormalization group 
evolutions of the soft function and the jet functions can combine to give local running. 
Taking a Fourier transform, 7(7/) = f dt exp{—ity)^{t) and F{y) = f dt exp(—ity)F(t) we 
have a simple multiplicative RGE 
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where y = y — iO. Note that the form of the position space anomalous dimension in Eq. ( |D6| ) 
simply follows from locality applied to the bi-local vacuum matrix element defining the jet 
function in position space. When translated to momentum space this directly implies the 
convolution structure shown in Eqs. ( |o7|Jo7| ). Integrating Eq. ( |D(j| ) from /i to /i by changing 
variables to a s with <iln/z = da s //3[a s ] gives the solution 
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(Note that in the main body of this paper where we consider LL and NLL accuracy we write 
to, K, and instead of uj, K, and K^, respectively.) Thus the position space solution is 
F(y,fi) =U(y,n,fi )F(y,no) with 
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The desired solution is the inverse transform, F(t,n) = 1/(2tt) fdy exp(ity) F(y,n) so 
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To simplify this result we use F(y,Ho) = J dt' exp(— it'y)F(t', Ho) and also the inverse trans- 
form 
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Doing the integrals over y and f we obtain the final result 
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where the evolution kernel is 



U(t- t',/i,/i ) 



(t-t>) 



(D13) 



This +-function is defined by Eq. ( |Cj] ) and K and u are determined at what ever order one 
desires from Eq. (|D8|) . In section |II [D|we carry out these integrals with NLL accuracy. The 
above derivation also suffices to solve Eqs. ([54]) and (|64|) to obtain Uh q and Uu m respectively. 
First we note that Jh q (Q, h) has the same form as ■yp(y,fi) with j = 2 and iye 7B — >■ 1/Q 2 , 
so the general solution is given by Eq. (|D9|) with the same substitutions, yielding the result 
in Eq. (|79|). For "iH m {Ql m i H) the cusp angle is fixed at m 2 /Q 2 , so the solution is given by 
Eq. (|D9D with r[a a ] -> and 7 [a s ] -> r Hm N ln(m 2 /Q 2 ) + 7Hm [a,]. This yields U Hm in 
Eq. 

For the cases with convolutions a few additional identities are useful. The evolution 
kernels obey 
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which states that it is equivalent to evolve through an intermediate scale, Ho — > Hi — > Hi or 
directly from Ho — > /J. To verify Eq. ( pi4[ ) one needs 
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(D15) 



where here u)\ = uj(hiHi)i u 2 = ^(hi->Ho)-i an d u ' = ^(Hi Ho) — %i + <^2- The first result in 
Eq. ( pi5|) is straightforward to derive using the Fourier transform. Another useful identity 



S3 



simplifies the convolution of two U's that have the same renormalization scales, but variables 
with different mass-dimension, and different anomalous dimension coefficients 

Jdr' U(Q'(r-r'),fJ l ,fJ, ;j',T',j',uji) U(r'-r", fj,, ^j, T, j, o5 2 ) 

Here the variables after the semicolon denote parameter dependence, and Q' simply denotes 
a variable with mass dimension j' — j. Also here U\ = u;(/i, /io; T'/f) and 0J2 = £>(/i, /io; T/j) 
are simply the u?s obtained from the other parameters, but this is not the case for uj' in the 
U on the RHS where uj' = uj\ + uj<i- The final useful identity is 

]im U(r-r',n,iio-J,T,j,u') = e R ^°' Tn) 5{r-r') , (D17) 



which is easy to derive from Eqs. ( pi3| ) and ( |CT|) . Using Eqs. ( pi5| ) and ( pi6| ) it is a 
straightforward exercise to verify the consistency equations directly in the integral form 
given in Eqs. ([71]) and (|77|). 



APPENDIX E: ANALYTIC RESULTS FOR G± IN THE PEAK & TAIL CROSS- 
SECTION 



In this appendix we show how the functions G± defined in Eq. ( |184j) can be determined 
analytically. As a first step we use Eqs. ( |142| ) and (|115| ) with k 3 = /ir and k 2 = firmj/Q to 
compute the £' integral to 0(a s ) [recall mj = mj(/ir)] 
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where z = sf /ir and 8± is given in Eq. ( |1 75| ) . The S'(z) term in Eq. ([El]) contains the residual 
mass correction 5mj for the jet mass scheme in Eq. ( 161 ) and the subtraction Si for the soft 
function from Eq. ( |175| ). Using the relations 
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the result for E± can be rewritten as 
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Note that for /xa ~ A and /ir ~ s ~ T ~ QA/m there are no large logs in this expression, and 
that the terms with 5m j and 5± are the same order in the power counting. Given Eq. (|E|). 
doing the second integral in the variable s" involving the Breit-Wigner function is simple 
since it just results in a shift of the invariant mass variable into the positive complex plane, 
as in Eq. fl37|): 
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(E5) 



For the final integration in the variable s' in Eq. (|184f ) we have to convolute the finite width 
version of the terms in Eq. (|E4j) with the evolution kernel Ub- The relevant computations 
read 
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with K 3 = K 3 (fi,fi ) and uii = Ct>i(/x, /Xo) given in Eqs. ( |147| ), and 
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Here H(u) is the harmonic number function and ty'(z) = d/dz\T'(z)/T(z)] is the derivative 
of the polygamma function. Thus the final result for the function G± is [rrij = mj(fi-p)] 
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APPENDIX F: CROSS-SECTION IN THE ULTRA- TAIL REGION 

In the text we presented results for the cross section in the peak and tail regions, where 
we assume Q ^> m ^> s and m ^> T. In this section we analyze the cross section for 
Q 3> m ~ s, so that s is far above the peak region. This corresponds to \M tj t — mj\ ~ m, 
where the top-antitop jet invariant mass double differential distribution can be described by 
the SCET factorization formula of Eq. ([13]). Writing this cross-section in an analogous form 
to Eq. (|180| ) [mj = mj(/i m )] we have 

d 2 a _ 4:M t M - t a rSCET( 
dM t dMf (mjT 

where 



^F sc ^(M t ,M { ,Q, mj ). (Fl) 



F sc ^{M u M- u Q,mj^ = (mjT t ) 2 H Q (Q,fi m )M(mj jf i m )uP Q (Q^ m ^ A ) 

/CO 
d£ + d£-G n {s t - Qt, Q, mj, T, fM A )G n (si - Q£~, Q, mj, T, // A ) S mod {£+, £~) . (F2) 
-oo 



Note that Eqs. ( |F1|) and (F2| ) are also appropriate for describing the massless limit mj — > 
and the stable limit T t — > 0. To obtain this result we manipulated the first form given in 
Eq. (|4"9D . Here, S mod is the hadronic model function given in Eq. ( |170| ), where we have 
suppressed its arguments a, b, A. The functions G n ^ n can be written as 



G nfi (s,Q,mj,T, fx A ) = J ds'ds"d£' Uj(s - s', /x A , /i r 



x J n ,n(s' ~ s" - Q£', mj, fi m ) 5 part (f , ha, <*i) , ,^2™ ■ ( F3 ) 
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Here, S' part is the modified partonic soft function of Eqs. (|115|) and ( |1 79|) . For consistency 
with our peak cross-section results we continue to use the jet-mass scheme, by taking 5mj 
for the Sm term in Eq. ( |107| ). Since it does not require any technical effort, we include 
in G n> n a constant width term for the top quark through the convolution involving the 
variable s" . We note, however, that away from the resonance region this width term leads to 
power-suppressed effects, and, moreover, does not provide a consistent description of the top 
quark decay. It is nevertheless convenient to introduce the width term for practical purposes 
because it allows for an easy numerical evaluation of the SCET factorization theorem for all 
values of M 4 f without running in to singularities for M t j close to the top quark mass. 

We can carry out an analytic calculation of the functions G n ^ defined in Eq. (|F3|) . The 
calculation divides itself into two parts, the terms singular for stf — ► which include the 
^-function and +- functions in Eq. (|106|) , and the non-singular ^-function term on the last 
line of Eq. (|106|) . The final result in the jet mass scheme reads [mj = mj(fi m )] 



G nt n (s, Q, mj, T t , fiA = Im g n>n (s, Q, mj, T t , /x m , /iA + 5Q nfi (s, Q, mj, T t , /i m , /z A 

■ ^nonsing / /-\ \ 
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For the singular terms, Q n ^ and <5{? n ,n, a computation can be carried out in close analogy 
to the computation of G± described in the App. [Ej Taking the singular terms in Eq. ( |106[ ) 
with k\ = fMrnin j, and the results from Eqs. ( |115| ) and ( |179| ) with k 2 = fi m m j/Q we find 
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Here, lo\ = u>i(fi\, /i m ) and K\ = Ki(fiA, Hm) are given by Eq. ( |105|) and the functions I n 
were denned in Eq. (|E8|) . The term 5Q nyfl arises from the jet mass definition (see Eq. ( |161| )) 



and from the subtraction that we carry out for the soft function model (see Eq. ( |1 75| ) ) . 

For the non-singular terms we note that the top-width effects represent 0(T t /m) power 
suppressed terms for any s <C Q 2 . For instance, when s ~ m? the top width appears in the 
combination s + imjT t , and when s ~ mT the entire non-singular term is 0(r t /m). Thus 
it is consistent to neglect the top quark width for the non-singular terms. Setting T t = in 
Eq. (F3) and carrying out the s' and s" integrals gives 
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APPENDIX G: ANALYTIC RESULTS FOR THE P T FUNCTION FOR THRUST 



In this appendix we derive an analytic result for Pt, the perturbative corrections appear- 
ing in the thrust cross-section in Eq. (|188|) . Starting from Eq. (|189|) for Pt with Eq. (|183[) 
for P the key is to simplify the integral over the product of G + G_, 
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where -B± ee (s, T t ) is simply a Breit-Wigner as shown in Eq. (p8|) . This Breit-Wigner appears 
due to the factorization of lifetime effects in section pi I Bj , and the perturbative corrections 
to the jet-functions are part of E± =0 . 

Here we rewrote the original integrations over symmetric variables (s s 's) and antisym- 
metric variables (s d s). Each of the integrations in square brackets can be performed, with 
the help of Eq. ( P16| ) and (|E1|) , to give the terms in the following result 
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Here the Breit-Wigner has a width of 2T t , and the function E^ =0 (s) is identical to E+ =0 (s) 
in Eq. (|E1| ) but with the replacements {a s — > 2a s ,5i — > 25i,5mj — > 25mj}. Finally the 
evolution kernel is 



U B (s, (i, Ho) 



n r(-2ti>i] 



(s - s') 1+2wi 



(G3) 



which is equivalent to the kernel Ub for the bHQET jet function but with the anomalous 
dimensions doubled. We can write the function E^r°(s) = Im[£!p =0 (s)] , where £^ =0 (s) 
has the same form as £+ =0 (s) in Eq. ( |E4j ) but with the replacements {a s — > 2a s ,5i — > 
2Si,8mj — > 25mj}. The result in Eq. ( |G2[ ) has a structure such that we can perform 
the last two integrations with the same techniques as in Appendix |E]. We obtain G^(s) = 
Iid.[Gt(s) + 5Gt(s)] where 
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and Gn(s) = G^(s, T t , //a, Hr) is simply G n (s, T t , /iA> A*r) from Eq. ([Erf ) but with u;i 
K 3 -> 2^3, and r t -> 2r t . 

Including the prefactor from Eq. (|189|) we have the final result 
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